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ANNIVERSARY OF THE SCIENTIST 


Yury Olegovich Chernyshey, engineer, teacher, scientist, is 85 


July 11, 2022, marks 85 years since birth of Yuri Olegovich Chernyshev, the Honored 
Scientist of the Russian Federation, Honorary Professor of DSTU, Professor of the 
Production Automation Department, Don State Technical University (DSTU), chief 
researcher of the R&D Center, Krasnodar Higher Military School named after Army 
General S. M. Shtemenko. He devoted 63 years of his life to scientific, pedagogical and 
scientific-organizational activities. 

Yuri O. Chernyshev graduated from Taganrog Radio Engineering Institute in 1959, 
defended his candidate's thesis in 1967, and his doctorate dissertation — in 1984. 

From 1959 to 1980, he worked in Rostov-on-Don. He was associate professor in Rostov 
Institute of Agricultural Engineering, as well as senior engineer — head of the group, 
Rostov Higher Military Command and Engineering School of Rocket Forces named after 
Artillery Chief Marshal M. I. Nedelin. From 1981 to 1983, he was Director of the North 
Caucasus Branch of the Russian Branch of the State Design and Technological Institute for the Mechanization of 
Accounting and Computing Works of the Central Statistical Bureau of the USSR. 

Since 1983, Yuri Olegovich has been Head of the Electrical Engineering and Automation Production Department, 
Rostov State Academy of Agricultural Machinery (RGASKHM). In 1988, he founded the Department of Applied 
Mathematics and Computer Science at this university and headed it. From joining RGASKHM to DSTU in 2009 to 
the present day, Yuri O. Chernyshev has been a professor at the Production Automation Department, Don State 


Technical University, and a chief researcher of the R&D Center, Krasnodar Higher Military School named after 
Army General S. M. Shtemenko. 

Yuri Olegovich Chernyshev is a leading scientist in the field of optimal design of production processes and computing 
structures, practical creation of computer-aided design systems, automatic process control systems. 

The scientist made a great contribution to the training of highly qualified scientific and pedagogical personnel. He 
enriched science with major scientific works and developments in the field of computer technology, theoretical and 
technical cybernetics. Under his leadership, about 80 fundamental state-budgeted and self-supporting R&D of national 
importance with a significant economic effect, confirmed by the relevant implementation acts, were carried out. 

Yuri Olegovich has created a successful scientific school in the field of non-traditional models and methods for solving 
optimization design problems. Nowadays, the school is still running successfully. He has trained 7 doctors and 35 
candidates of sciences. 

His fruitful scientific activity is presented by more than 500 research papers, including 17 monographs, 22 patents, 18 
programs in the funds of algorithms and programs. 

Professor Yu. O. Chernyshev was recognized by the domestic and foreign scientific community. From 1976 to 2016, he 
was a member of three dissertation councils for the defense of doctoral and candidate dissertations. From 1995 to 2000, he 
was a reviewer of the international journal “Mathematical Reviews” (USA, University of Michigan); a member of the 
editorial board of the journal “Advanced Engineering Research”, a member of the American Mathematical Society. 

In 2001, he received a Letter of Acknowledgement from Governor of the Rostov region “For a great personal 
contribution to the development of Don science and vocational education”. 

In 1996, by Decree of President of the Russian Federation, Yuri Olegovich was awarded the honorary title “Honored 
Scientist of the Russian Federation”, in 2003 he was awarded the Order of Friendship, in 2008 — the Vernadsky Gold 
Medal, and in 2020 — the L. V. Krasnichenko medal. In 2011, Yu. O. Chernyshev was given the title of “Honorary 
Professor of DSTU”. 
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Abstract 

Introduction. The article studies the problem of investigation of coupled nonstationary thermoelectroelastic fields in 


piezoceramic structures. The main approaches related to the construction of a general solution to the initial non-self- 
adjoint equations describing the process under consideration are briefly outlined. The work aims at constructing a new 
closed solution to the axisymmetric thermoelectroelasticity problem for a long piezoceramic cylinder. 

Materials and Methods. A \ong hollow cylinder whose electrodated surfaces were connected to a measuring device 
with large input resistance was considered. On the cylindrical surfaces of the plate, a time-varying temperature was 
given. The hyperbolic theory of Lord—Shulman thermo-electro-elasticity was used. The closed solution is constructed 
using a generalized method of finite integral transformations. 

Results. The developed calculation algorithm makes it possible to determine the stress—strain state of the cylinder, its 
temperature, and electric fields. In addition, it becomes possible to investigate the coupling of fields in a piezoceramic 
cylinder, as well as to analyze the effect of relaxation of the heat flow on the fields under consideration. 

Discussion and Conclusion. The use of assumptions about the equality of the components of the temperature stress 
tensor and the absence of temperature effect on the electric field allowed us to formulate a self-adjoint initial system of 


equations and construct a closed solution. 


Keywords: thermoelectroelasticity, hyperbolic theory, nonstationary coupled problem, long piezoceramic cylinder, 


finite integral transformations. 
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Introduction. Recently, different-purpose technical devices made of piezoceramic material have become 
widespread. Here, devices whose operation is based on the effect of the coupling of elastic, electric and temperature 
fields hold a specific place [1]. Various theories of thermoelectroelasticity were developed to describe their work taking 
into account the coupling of fields [2-4]. At the same time, for a better description and evaluation of non-stationary 
processes in structures, it was required to construct analytical solutions. However, the mathematical formulation of the 
problems under consideration included a system of non-self-adjoint partial differential equations, whose integration was 
difficult to treat mathematically. 

To solve this problem, as a rule, equations are investigated in an uncoupled form [5, 6], infinitely long bodies are 


analyzed [7—11], or thermoelectroelasticity problems are considered in a quasi-static formulation [12, 13]. 
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In this paper, we consider a coupled dynamic thermoelectroelasticity problem for an infinitely long hollow 
piezoceramic cylinder. As a result of the transformation of the initial calculated ratios, it was possible to form a self- 
adjoint system of equations, the integration of which was carried out by the method of incomplete separation of 
variables in the form of a generalized finite integral transformation [13]. 

Materials and Methods. Let a hollow, long, loose in the radial plane, piezoceramic cylinder occupy area © in the 


cylindrical coordinate system (r,. 0, z) : {a <r,<b,0<S0<21,-10<z< oo} . On the cylindrical surfaces, the 
temperature is given in the form of the following nonstationary functions (boundary conditions of the Ist kind) — 
@:(t.) (4, =a), 5(t.) A =); The internal electrodated surface is grounded, and the external one is connected to a 
measuring device with a large input resistance (electric idle mode. 

The mathematical formulation of the axisymmetric problem under consideration in dimensionless form includes 
differential equations of motion, electrostatics, and thermal balance based on the hyperbolic Lord—Shulman theory, as 


well as the boundary conditions [2, 7, 15]: 


2 
Gg ay yg a youn! 29, (1) 
Or r- Or r Or “r Ot 
Vv bogges ee aves: 
r Or r Or 
oO [A Oo 
—-| — a,O+a,VU )=0; 
or (4 po VU) 
OCU U & 
r= Rl +4) — + = -O=0 Oj. 4 =O}, O4 = Ors (2) 
a0) OU 
sp =0., + +a,—+a,0 =0; 
tyr ( Be B¥ 5 6 |. 
ou -0, 0, =e. (3) 
Ot Ot 
' Gia ah Wipe. 3 : ea Peat 1G C 
where {U,r,R' =}U",r,,as/b ,d=—*-@ , {0,0,,0,'=—30"",o' -T,, 0, -T,}, {t, B\ = 3 a-—, 
{U,r,R} ={ \ args ere | ~ Ty; ~Th} {t,B} =~ oe. 
a,=—' 4,-™, a,= a a, = 3 a ee eee bY sT ‘ie 
53 33 Cy3€33 C3333 €33733 AN\ p AYC,3p Cy, 


U'(r,,t.),0 (7,t.),O' (%,t.) — accordingly, the radial component of the displacement vector, the electric field 


potential, and the temperature increment in dimensional form; ( ©" (7,,t,)=T7(7,,t.)-T)(%), — current temperature 


as? P> mss €33— elastic moduli, density, piezo module, and 


t 


and temperature of the initial state of the body; C 
permittivity coefficient of electroelastic anisotropic material; (m,s = 3 V11>7%33— components of the temperature 


stress tensor (7,, =C,,@,,7%3; =C,,a@,); A,k,a,— coefficients of thermal conductivity, volumetric heat capacity, and 
linear thermal expansion of the material; g,— component of the pyroelectric coefficient tensor, £,,,— relaxation time; 


: ie 0 1 
©,— the rate of temperature change known at the initial moment; V = ap +-. 
r + 


In the case of grounding of the inner surface of a piezoceramic element, electrical voltage V(t.) is determined by 
the potential on its outer surface: 

V(t,)=¢(Lt.). (4) 

When constructing a general solution at the first stage of the study, the radial component of the electric field 


intensity vector is determined as a result of integrating the electrostatics equation: 
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U D, 
+d,—+a,0+—, (5) 
r r 


where D, — integration constant. 
Substituting (5) into (1)—(3) allows us to formulate a new problem with respect to functions U(r,t),@(r,t) .In 


this case, the condition of the absence of a radial component of the electric field induction vector on the outer 


cylindrical surface of the element (the last equality (2)) is fulfilled at D, =0; and the condition of grounding of 
the inner surface (Gr = 0) is satisfied as a result of integration (5). 


At the next stage of the solution, the inhomogeneous boundary conditions (2) are reduced to a form that allows 


further use of the procedure for incomplete separation of variables by the method of finite integral transformations. To 


do this, new functions u(r,t), N(r,t) related to U(r,t) ,O(r,t) are introduced: 
U(r,t)=A,(r,t)+u(r,t), O(7,t) =A, (r,t) +N(r,1), (6) 
where H,(r,t)=f,(r)A(LA+ 46 (r)A(R1)+ A (NO (+h (r)a, (0), 
H,(r.t)= f.(r)@(1)+ fo(")@, (1), 
f,(r)...fe(7) — twice differentiable function, A(r,t)_,, =(1+a,—a5—a,)U(r,1)/r- 


r=R,1 
Substituting (6) into the calculated ratios (1)-(3) with respect to functions U (r,t),O(7,1) when the following 
conditions are met: 


(1+a,)VH, +(a,—-1)H, = A(r,t) Hy,-2 =, Hy,. =, (7) 


r=R,1? 


allows us to get a new boundary value problem with respect to function u(r,t) ,N (r,t) : 


Ou u ON N 1 Gu 


“OF OR os Chal ee “ 
vE (Sop SZ \anravu)=F 
r=R,1Vu=0,N=0; (9) 
t=0u=-H,(r,0) ,.N=-H,(r,0), (10) 
6u__ OH, (r,t) ON _6 GH, (t,t) | 
Ot ot °or ° or” 


where F, = v= + a + pe lem. +ava) 
r 


H H HH H 1 °H 
F = yey ayaa 2 aI 
Or r Or ae 3 (1+a,) ot 
+ -1 
p=" a,a; poet dha a, 
+a, l+a, l+a, 


It should be noted here that A(r,t) is a function of the displacements of the cylindrical surfaces of the cylinder. 


Initially, A(r,t) is equated to zero with its subsequent determination and refinement H,, F,, F,. 


Further transformations of the calculated ratios (8)—(10) are associated with the use of the following assumptions: 


b, =1,b, =0, and the introduction of a thermoelastic potential 


N=VB. (11) 
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Condition b,=1 can be accepted without a large error, since for piezoceramic materials b, =0.94+0.98, and 
dependence b, =0 is fulfilled in the case of equality of the components of the temperature stress tensor (7,, = 73, ) and 
the absence of temperature influence on the electric field ( g, =0). 


As a result, the following task is formed regarding u(r,t),B(r,t) ; 


1 ou 


o o 
—Vu+t+b,—VB =F, 12 
or ut * Or (1+a,) ot” a 
o 0 e 
reac [S66 Z |laaeau) =F 
rRivesve=0: (13) 
t=0 u=-H,(r,0) ,VB=-H,(r,0), (14) 
Ou _ antag 2 yp 6, 2), 
ot ot ot Ot 
OH 0 a 
where F, = ae {4 pF Nan +a,H,),VH, = f;(r)@,(t)+ f(r) @ (t). 


The initial boundary value problem (12)—(14) is solved using the structural algorithm of the generalized finite 


integral transformation (FIT) [14]. At the same time, it is possible to use single-component unknown transformation 


kernel K (Ay r) for this task: 


{G,(/,,1),G,(A,,)} = [{u.), BU. )}KG, rar, (15) 


{u(r,t),B(r,t)} = iG (A,.t),G)(A,,.01K (Ar) Kl. (16) 


=1 


1 
|K, |’ = [KG,, r)’rdr ; 
R 


where A, — eigenvalues forming a countable set. 


As a result of using the FIT algorithm [14], we obtain problems with respect to the transformation kernel K (A, r) ; 


°K(A, IK (A, 
g Cee ( wD ala = ]KA.N=0, (17) 
dr r dr r 
r=R1VK(A,r)=0, (18) 
and transform G, (/,,t),G,(/,,1): 
ae 1 dG, 
A?G,, + —G,, ad eee (19) 


2 d d* . 
A; G,,; ne ae (a,G,, +4,G,,) = Foy > 


dG,. : dG,, : 
t=0G, =Gip, ap = Gin Gr; = Gyo» a = Gyi9; (20) 


1 1 
where {Fiy,Fy,}=[{F.F}K(A,r)rdr, {Gos Goro} =—-[{H, (1.0), H, (7,0)} KA, rar, 
R R 


[t= ff (e.- A) |aon 
R ° |t=0 
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The general solution to problem (17), (18) has the form: 
K(A 1) =% (A) A (Ar) —Jo (4 )¥ (Ar). (21) 
Here, eigenvalues A, are determined using the following transcendental equation: 
Y (A) Jo (A,R)-J, (4,)% (AR) =0. 
The system of differential equations (19) is reduced to the following resolving equation of the 4th order with respect 
to G(4,,0): 


d* d° d* 
+b 
= “dt dt? 


+s; + Dg, “ +b, |G, =F, (22) 
t 


(1+a,)( 2 o 
where F,, =-—— J lh +(i+a,) )Fu | z [S16 S hn, 
2 1+ 
vetinee[ienet } =4 Goer =a! a4) . 
B a, a,B B a, a,p 


Since the characteristic equation corresponding to (22), 
k* +b, +b k? +bk+b,, =0, 


is valid, then it, from the condition of the oscillating solution for G,(/,,f), has two real roots (k,,,k,,;) and two 


li? 
complex-conjugate roots 
(k,,=o0+i0, k,, =o-i@). 
In this case, the general solution to equation (22) has the form: 


G,(/,,t) = D,, exp(k,t) + D,, exp(k,,t) + D,, exp(k,t) + Ds, exp(k,,t)+ (23) 


wl Fa (r){exp[k, (t- t)| —exp| k,, (t —r) |} dt by, (z)exp| o(t 7) {Pro sin(wt—ar)—b,,,cos(at+ar)}dr, 


-1 


where b,, = {i —k,; i, 7 o) +o" } 


bh, =[@ (03, +03,)] 0-240 


10i 


= kk; -(k; tk; )oto’-o’, 


k,, -k 


2i 


b,, = 0(20-k, -k,,),@= 3§ Kai} 


Function G,(4,,f) is determined from the first equation of system (19). Substitution of the obtained 
expressions for the transforms under the boundary conditions (20) makes it possible to determine the integration 
constants D,,...D,;... (19). 


Substitution of G,(A,,t),G,(/,,f) into (16), (11), (6) allows us to get the final expressions for functions U (r,t) ‘ 


O(r,t): 


—2 


U(rnt)=H,(r.t)t 0G AK (ANIKI, (24) 


i=l 


O(r,t)= + La, VK (A,r) ||K, |. 


At the final stage of the study, functions H, (752), H, (7,1) are determined through solving the following differential 


equations: 


Fg epg SG, (25) 
or r or 


which makes it possible to significantly simplify the right parts (F, F,) of the calculated ratios (8). 


Mechanics 


85 


http://vestnik-donstu.ru 


86 


Advanced Engineering Research 2022. V. 22, no. 2. P. 81-90. ISSN 2687-1653 


Substitution of expressions for H,,H,in (25) enables to form systems of equations with respect to functions 


J, (r)...£,(r) , that are determined when the conditions are satisfied (7). 
The potential of the electric field of a piezoceramic cylinder is determined from integrating equality (5) and 


satisfying the next-to-last boundary condition (2): 


male OH, (r,t) +a, H, (r,t) +a ( rt) fr (26) 
or r 


* + De (t) 


K; 


+>GA,,0)B(A)||K | +4. GA,.0B(A)|K, 
i=l i=l 


where D,(t) = J c atid ep (44) +a,H, (-1) fre Sa (A,.0B, (AIK, + 
r i=l 


4G NBLANIK | 1m (4)= a OK (4,1) +a, Se, 


= or r 


B,(A)=JVK(A,r)er. 


The obtained calculated relations (24), (26) satisfy differential equations (1) and boundary conditions (2), (3), i.e., 
they are a closed solution to the problem under consideration. 


Research Results. As an example, we considered a radially polarized piezoceramic cylinder (b =0.02m, R=0.8 ) 


of PZT—4 composition, having the following physical characteristics [10]: = 7500 kg/m* , 


{Ci1sCyssC,s} = (13.9; 11.5; 7.43} x10" Nim?, {e,,,€,}={-5.2; 15.1} Chm’, 
{Hiss Yeat = (4.6; 3.9} 10° H/(m* °C), €,, =5.62x10°’ ®/m, g, =2x10*Ka/(m’ °C), 


k =3x10° J m?* °C), A=1.6 W/(m °C), B,,, =5x10°c. 
A temperature load acts on the inner surface (7, =a ) of the piezoceramic cylinder: 


a, (Aly ae in( Et | -1,)+H(t, — tr.) Oo; (t,) =0, 


where H(?) —Heaviside step function (H(7)=1 at 720, H(f)=0 at 7 <0), 7, =Te —Tys Trhagsl nox — 
maximum value of the external temperature effect and the corresponding time in dimensional form 
Ge, 3373 K (100°C), 7, = 203-K (20°C), F.=1 ©), 

Figure 1 shows graphs of changes in functions ©* (a), U (24), (r,t) along radial coordinate r at various points 


in time t. The numbers 1—3 respectively indicate the results obtained at the following time values: t =1,,,.54t ax >! Snax 


Analysis of the calculation results allows us to draw the following conclusions: 

— sufficiently large value of the coefficient of linear thermal expansion a@, of the piezoceramic material causes rapid 
heating of the cylinder; 

— radial displacements on the inner cylindrical surface (r=R) at the first stage of the study (t=1,,,,) take the 
greatest values, followed by a decrease over time. The reverse pattern is observed with respect to the displacements at 


r=1; 


ry 
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@' (7,1), °C 


c) 
Fig. 1. Graphs of ©* (7,1), U(r,t),A(7,t) changes along radial coordinate r 


at various points in time t (1—t,4.52—Afinax >3—15tingx ! d— O" (r,t)+r 3 b— U(r,t)=r;c— O(r,t)+r 


The degree of coupling of thermoelectroelastic fields is most conveniently analyzed using the coefficient 


2 
by, = z( a,t+ 2 from equality (21). Here, a, determines the coupling of electroelastic fields, and “ _ the 
a, a, 
effect of the rate of change in the volume of the body on its temperature field. 
Figure 2 shows displacement graph U (L t) in time f taking into account (solid line) and without taking into 
account (dotted line) the induced electric field. 
It should be noted that the preliminary polarization of piezoceramics causes the formation of a more “rigid” material 


(a, = 0.353) and, accordingly, a decrease in displacements during deformation of the cylinder. 
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The coupling of temperature and electroelastic fields in a piezoceramic cylinder can be neglected due to the small 


value 


%s -1.8x107 «1. 
a, 


U(r,t) «10° 


0 1x10° 2x10° 3x10° 


Fig. 2. Graph of change U (1, t) in time ¢ 


(solid line — including the polarization; dotted line — without including the electric field) 


Figure 3 shows graphs of changes in electrical voltage V(t) over time, taking into account (solid line) and without 


taking into account (dotted line, £,, = 0) the heat flow relaxation. 


V(t) x 10° 


Fig. 3. Graph of change V(t) over time t(t,,, =0.001 s,t,,. =196) 


The calculation results show that for the problem under consideration, the refined hyperbolic Lord-Shulman theory 


a) >5.56x10° K/’and at 
dt. s 


|*=0 


should be used at a high rate of change in the temperature load ( t.,, > 0.001 s, 


lower speeds — the classical theory of thermoelectroelasticity ( £,,, =0). 

Discussion and Conclusion. The constructed new closed solution to the coupled dynamic problem with satisfaction 
of the boundary conditions of thermal conductivity of the Ist kind made it possible to determine all the components of 
thermoelectroelastic fields in a long piezoceramic cylinder. The advantage of the presented calculation algorithm is that 
there is no need to approximate the temperature function when studying the equation of motion, in contrast to the 
uncoupled formulation of the problem. At the same time, the effect of the rate of change in the volume of a 


piezoceramic body on its temperature field can be neglected. 
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Abstract 

Introduction. The paper presents a brief kinematic analysis, as well as the application of D'Alembert's principle to 
finding the relationship between the force parameters in the transmission mechanism of a robotic manipulator 
constructed from a twisted arm chain. The use of this transfer mechanism can enhance the life of the arm actuator, the 
accuracy of its positioning, and increase the workload compared to the flexible linkage actuators (twisted strings). The 
work aimed at obtaining dependences between the displacements of circuit elements, as well as their force parameters 
required to monitor the control system operation of these devices. 

Materials and Methods. In the course of solving the problem, an elementary segment (element) of the transmission 
chain was considered. To find the relationship between the loads in the element, the virtual displacement principle was 
used. When finding kinematic connections between displacements, a brief geometric analysis of the transmission chain 
element was carried out. To conduct a comparative analysis of the analytical dependences obtained, a simulation 
technique implemented on a graphical link model using the NX software package was applied. 

Results. In the course of the study, we obtained dependences for determining the magnitude of the moment developed 
on the input link, depending on the external workload and its rotation angle, as well as for defining the linear 
displacement of the output link. A simulation model of the actuator was constructed, which can be applied in the 
dynamic study of the actuator mechanism, taking into account the inertia of the links. 

Discussion and Conclusions. From the obtained analytical dependences, we determined the value of the angle of 
rotation of the input link of the mechanism element, at which the maximum torque value for a fixed workload on the 
output element was achieved, as well as the maximum linear displacement of the output link. The calculated values 
were in good agreement with similar values obtained from the results of the simulation experiment, which gave us the 
possibility of using analytical dependences in the formation of a robot control system. In addition, these dependences 


made it possible to provide the selection of actuators with the required force indicators. 


Keywords: manipulators, flexible linkage, driven elementary segment, rigid arms, virtual displacement principle, 


torque, simulation experiment. 
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Introduction. Recently, devices based on twisting cables or strings (threads) have begun to gain popularity. These 
inventions are universal in use and are in demand in various areas of robotics, specifically, in the manufacture of 
manipulators, grippers and exoskeletons’ [1-4]. Drive structures based on twisting flexible elements have a number of 
advantages. These include the size of the device, ease of assembly and installation, versatility. As an example, a 
mechatronic twisted actuator, which is an electromechanical device with numerical control, is able to independently 
monitor the state of the working body based on data obtained from the shaft position sensor and the force sensor. 

In the presented work, the authors carried out a kinematic and power analysis of the transmission mechanism of the 
actuator, built on the basis of a twisted arm chain. This design is intended to replace actuators using flexible 
transmission links, namely, twisting cables and strings. 

The experimental device shown in Figure | consists of an electric motor and several strings oriented along the motor 
rotation axis. Strings connect the motor drive shaft to the simulated load generated by another electric actuator. Load 


control (linear actuator) enables to consider this system as a spring damper with adjustable parameters. 


Linear 
actuator 


Force sensor Rotary actuator Position 
sensor 


Fig. 1. Experimental setup [4] 


Thanks to the string-twisted actuator, it was possible to place all the drives and motors in the prosthetic arm, given 
that all the main degrees of mobility of the arm were provided [5—8]. At the same time, the weakest point of the system 
is a flexible cable or string. Wear and stretching of the twisting string are the key disadvantages of this type of actuators. 

An alternative solution to this problem may be the replacement of twisted strings or cables with a chain link system. 
Let us replace a fragment of twisted strings or cables with an element (segment) of rigid levers (rods) that pivotally 
connect two circles with radius R and centers O; and O, (Fig. 2 a). Input link / is supported by thrust bearing 6 and has 
the possibility of rotational movement relative to center O;, located on Zp axis. In this case, lower circle 4 (the output 
link) can make a translational movement along Zp, axis due to the presence of guide 5. As a result of rotation of the 
upper circle, levers 3 and 4 perform spherical motion relative to link /, thereby changing distance / and moving 


workload P (Fig. 2 b, c). 


: Gaponov II, Nedelchev SI. Mechatronic drive on twisted threads. RF Patent no. 2020144036, IPC F 16 H 35/00; 2021 (In Russ.) 
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Hel 


P 
a) b) c) 


Fig. 2. Diagram of the actuator chain segment: a — in initial position; 


b — in intermediate position; c — in limit position 


For the practical implementation of this design solution, it is required to have a control system based on a 
mathematical model of the actuator, which reflects the kinematic and power connections between the elements of the 
actuator chain. 

This work aimed at finding kinematic and power connections in an elementary segment of the actuator chain. 

Materials and Methods. We propose to find the solution to this problem using the virtual displacement 
principle [9], which enables to exclude internal forces in the links of the actuator element from consideration. Thus, the 
torque applied to the input link remains under consideration, as well as workload P applied to the output link. 
Considering the scheme shown in Fig. 4 b, we take rotation angle g as the generalized coordinate for the input link, and 
linear displacement z — for the output link [10]. 

In accordance with the virtual displacement principle [11], we make up the equation of works: 

Mo6¢-P6z=0, (1) 


where 6g — rotational virtual displacement of the input link, 6z — linear virtual displacement of the output link. 
Distance H between links / and 2 is related to the linear displacement z of link 2, taking into account the length of 


the lever AB = /, dependence H =/—z. Then, taking into account the scheme shown in Figure 4 b, we can write: 


EF -(2rsin(£)) =(I-z). 


Value z, depending on the rotation angle q, is defined as: 


oi i -[2Rsin( $) . (2) 


According to [11], the relationship between virtual displacements can be established as: 


d 
OZ= fg) 1) : 
dg 
d 
Then, taking derivative a from (2) and performing the transformations, we get the dependence between the 
? 
virtual displacements 
R’ sing 


Oz= og. 


f-(xo(§) 


Substituting dz into the equation of works (1), we obtain: 
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R’ sing 


Zeno) 


whence, through reducing by 6g, we express moment M as rotation angle function g in the form: 


PR’ sing . (3) 


f-orse(§) 


To evaluate the obtained analytical dependences of moment (3) developed at the input link and the linear 


Mo¢-P 


56 =0, 


M= 


displacement of the output link from the angle of rotation (2), a simulation experiment was carried out [12, 13], 
implemented in the NX Nastran software environment [14]. For this purpose, a three-dimensional solid-state model of a 
section of the actuator chain (Fig. 3 a), consisting of two sequentially connected elementary segments of the mechanism 
under study, was constructed. One of the segments of the actuator (Fig. 3 b) was used to perform the calculation with 
the following parameters: R = 12 mm, / = 54 mm. When constructing the model, spherical joints were presented in the 
form of toroidal joints, which made it possible to perform chain elements in the form of curved rings connected by rigid 
levers (Fig. 3). This design allows the levers to perform spherical movement. Since this connection provides the levers 


with the required degree of freedom, it can be considered corresponding to a spherical joint (Fig. 5 c). 


L002 1003 1001 LOO L002 


(003 loot L001. | 


: >t! ~ 
L007 L005 L004 aft 4 tH 
NG, zx . cb 
1002) 
Ni ay 
a) b) @) 


Fig. 3. Three-dimensional solid models: a — chain section of two segments; 


b — chain segment; c — projections of a complexly curved ring 


In the developed three-dimensional model, kinematic connections were assigned, the centers of mass of the links 
were determined, and workload P = 100 N was set, after which the motion of the elementary segment of the actuator 
chain was simulated in NX Nastran [14, 15]. The parameters of the segment model are given in Table 1. 


Table 1 
Simulation Options 
Moment of inertia about Original coordinate of center 
Connection each of the axes, of mass, mm 
nuEADEY in kg: mm? Mass, kg Name 
simulation 1X Ty i X Y Z 
LOO1 0 1.63 209 ; 
L004 5.51 4.18 | 4.18 0.04 0 1.63 139.5 nee 
complex shape 
L007 0 1.63 69.8 
L002 —0.03 —12.9 173.4 
acel 0.03 12.9 | 173.4 0.02 — se = Lever 
L005 12.9 1.5 103.7 
L006 —12.9 1.5 103.7 
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Research Results. In accordance with the obtained analytical dependences for the amount of linear displacement z, 
as well as torque M, a calculation was performed in the MathCAD environment, whose results are shown in Figure 4. 
The calculation was based on the following initial data: workload P = 100 N, radius of the input and output link R = 12 


mm, lever length R = 12 mm, range of rotation angle g of the input link within 0—-180°. 


Zz, mm M,N-mm 
5.60 281.3 
ee 
4.70 —— 234.4 
ral 

3.75 va 187.5 

# 
2.82 Z 140.7 i 
1.90 / \ 

- 93.8 
0.94 J 46.9 iN 
/ 
0 30 60 90 120 150 180 @° 0 30 60 90 120 150 180 @,° 
a) b) 


Fig. 4. Dependences of the studied values on the angle of rotation: 


a — linear displacement of the output link; b — torque M at the input link 


From the presented graphical dependences, it could be established that the maximum displacement amount was 
achieved at a rotation angle approaching 180 ° and is 5.6 mm, and the torque data peaked at 281.3 N-mm at rotation 
angle g = 93°. 

As a result of simulation modeling, a displacement schedule of output link L004 was obtained depending on the 
rotation of the actuator link LOO1 (Fig. 5). The movement of the elementary segment of the lever mechanism occurred 
until the right and left levers touched in the extreme position. Since the levers had volume, a complete rotation of the 
output link was impossible due to the side surface contact, which is also shown on the graph (Fig. 5). To achieve a 
complete rotation, it is required to change the structure of the lever to a more complex shape. Figure 6 shows the 


dependence of the torque when turning the chain segment of the calculated model by 180°. 


Zz, mm 


144.9 


143.8 
142.6 
141.6» 
140.5 


139.5 


0 60 120 180 @,° 


Fig. 5. Travel schedule of the output link in the relative coordinate system 
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164.90 


82.90 


0 60 120 180 9,° 


Fig. 6. Graph of torque change during the simulation experiment 
The values of the studied parameters obtained through simulation (simulation experiment) differed from the 
values obtained theoretically by no more than 3 %. 

Discussion and Conclusions. Within the framework of the presented work, kinematic and power analyses of the 
elementary segment of the actuator chain of the manipulator were carried out. The obtained analytical dependences 
made it possible to determine the maximum amount of linear displacement and the developed torque in the actuator 
segment. The results of the simulation experiment implemented in NX Nastran environment and the subsequent 
comparative analysis showed a good correspondence of the obtained analytical dependences with respect to the three- 
dimensional dynamic model of the object. This allowed us to conclude that the obtained analytical dependences could 
be used in the development of a control system for the actuators of the above-mentioned manipulator devices. The 
conducted studies showed that the chain of segments of the lever mechanism in question had good transmission 
properties and, with sufficient length, could replace the transmission based on twisting the string. However, the device 
of the studied transmission is more complex. Despite this, the transmission based on twisting a chain of segments has 
the advantage of being able to provide a greater number of operating cycles of the manipulator device while 


maintaining the required positioning accuracy and repeatability due to less deformations of the chain elements. 
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Abstract 
Introduction. The article elucidates increasing the efficiency of turn-milling of powdered metal high-speed steel 


products. Turn-milling can be used as an alternative to the traditional turning method. The article describes advantages 
of the turn-milling method. A review of studies devoted to improving the surface quality of parts when turning by 
milling is given. The work aims at determining the effect of cutting modes on the surface roughness by the orthogonal 
turning method through milling powdered high-speed steel with a monolithic cutter. 

Materials and Methods. Statistical analysis methods based on the creation of a mathematical model for predicting 
microgeometric deviations of the treated surface were used. An experimental research method was applied to verify the 
adequacy of the mathematical model. The experiment was planned according to the non-composite design proposed by 
Box and Behnken. The experiment was carried out on a turning machining center with a driving tool. Powdered high- 
speed steel BOHLER S390 MICROCLEAN was used as sample material for the experiment. A monolithic carbide 
milling cutter served as a cutting tool. During the experiment, the cutting speed, milling width, and feed per tooth 
varied. The roughness of the treated surface was measured by a contact profilometer. 

Results. A mathematical model of the formation of surface roughness depending on the processing modes was 
developed. During the experiments, the effect of cutting speed, tool feed, and radial cutting depth on the roughness of 
the treated surface was determined. It was established that the dependence of roughness on feed had a linear character 
over the entire investigated range of cutting modes. In turn, the dependence of roughness on the cutting speed and 
cutting width had a parabolic character. The results obtained allowed us to achieve the roughness of the treated surface 
Ra = 1.85 without reducing the processing performance. 

Discussion and Conclusions. The developed mathematical model reflects the impact of cutting modes on the 
surface roughness when turning high-speed steel with a monolithic cutter. The results of the conducted research 
can be used to determine the optimal cutting modes that provide a given surface quality in the manufacture of real 
parts under the production conditions. It is recommended to continue the research with the control of additional 
output parameters, such as temperature and vibration. Reducing the effect of regenerative self-oscillations on the 
roughness of the treated surface can be reached through assigning the cutting modes based on the results of a 


modal analysis of the process system. 


Keywords: turn-milling, carbide cutter, roughness, high-speed steel, cutting simulation. 
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Introduction. An alternative to the classical cutting process is turn-milling. The first attempts to describe this type 
of processing and its practical application took place in the first half of the twentieth century. In 1948, G. I. Granovsky 
described the kinematics of processing a cylindrical rotating part using a milling cutter [1]. In 1964, A. O. Ehtin 
described various methods of turn-milling, including milling with a tool with an axis parallel to the axis of the part, and 
milling with a cylindrical cutter with circular and tangential feeds [2]. Schultz and Spur are considered to be the 
founders of the practical application of this method. They divided turn-milling operations into two groups: orthogonal 
and tangential [3]. A significant contribution to the study of this technique was made by U. Karaguzel, M. Bakkal, and 
E. Budak, who studied the mechanics and thermal aspects of turn-milling. They also developed models to test the 
geometry of the process, kinematics and mechanics, and the quality of machined parts in orthogonal and tangential turn- 
milling operations [4]. 

When implementing the turn-milling operation, a challenging issue is to obtain high surface quality, first of all, to 
reduce roughness. The analysis of publications allowed us to identify the main directions in the study of ways to 
achieve high surface quality: 

— study of the effect of eccentricity [5]; 

— determination of the profile of the cutting tool with multi-faceted indexable inserts with round and rounded 
cutting edges [6]; 

— study of the effect of temperature in the cutting zone [7]; 

— prediction of the wear resistance of the cutting edge when using the turn-milling operation [8]. 

It should be noted that in the above works, traditional structural steels were used as the materials under study. 
Investigations on high-speed steel, including those obtained by powder metallurgy, are not common. This limits the 
application of the advantages of turn-milling for the manufacture of rotating cutting tools [9]. This direction of 
machining can be attributed to the area, where turn-milling can give the greatest economic effect. 

The surface quality is affected by the cutting tool used. In particular, the application of end mills with replaceable 
polyhedral plates does not always provide the required surface roughness [10]. It is important to note that the considered 
experimental part of the research was carried out on multi-purpose machines. However, turn-milling can also be used on 
turning machining centers with driven tools [11]. 

This work aimed at determining the effect of cutting modes on the surface roughness in the orthogonal method of 
turn-milling powdered high-speed steel with a solid cutter. 

Materials and Methods. The experiment was carried out on a turning machining center with DMG NEF 400 drive 


tool (Fig. 1). 
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Fig. 1. Treatment zone during the experiment 

Powdered high-speed steel BOHLER $390 MICROCLEAN’ was used as the test material. The workpiece was a 
cylinder measuring 88x300 mm that was fixed in a three-jaw hydraulic chuck with preloaded center mounted in the 
tailstock of the machine (Fig. 2). Monolithic carbide milling cutter YG-1 GM999 12(R3)x1232x75, z =5 was used as 
a cutting tool. The treatment was carried out without the use of coolant. The roughness of the treated surface was 
measured using a stationary Taylor&Hobson Form Talysurf 200 profilometer. 

Chemical composition of BOHLER $390 MICROCLEAN steel: 

1.64 % C; 0.60 % Si; 0.30 % Mn; 4.80 % Cr; 2.00 % Mo; 4.80 % V; 10.40 % W; 8.00 % Co. 

Processing characteristics of BOHLER S390 MICROCLEAN steel, presented in a scoring system of 5 points 
corresponding to the maximum values of the properties of this group of steels: 

— compressive strength — 4; 

—  grindability — 3; 

— hardenability — 4; 


— viscosity — 4; 


— wear resistance — 4. 


'M390-BOHLER Russia. www.bohlernn.ru URL: https://www.bohlernn.ru/ru/products/m390/ (accessed: 14.12.2021). 
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Fig. 2. Scheme of orthogonal processing with YG-1 GMG19910 monolithic milling cutter 


The mathematical model of the deviation of the surface roughness can be presented by the equation [11]: 
InRa = InD + c, «Inf, + Cy, Inf? + cy: InV, +c. : InV2 +c - Ina, + C33: Ina?, 
where f, — feed, mm/tooth; V.— cutting speed, m/min; a, — cutting width, mm. 
Rename this expression as follows: 
V1 = Go Xo + G1 %1 + 911 XP +. G2 Xz + G22 XZ +G3°X3 + 933° XG, 
where y, = In Ra; x» — dummy variable; x), x2, x3 — coded factor values; go, gi, 92 g3 — constant coefficients. 
Encoding of variables is carried out according to the following relations: 
2dnf, +n fi nav) 
= — +1 
In Teint —In femin 
a 20dnY + Inv, 4: 
In Vespa —In Vemin 
2(In a, + Inde, ,) 
IN Ge nice — Genin 
— upper levels of factors; f,, 54» Vemin’ Lemin —~ lower levels of factors; f,,,Vc,,@e, — main 


1 


Xa +1, 


where Teaagh Viaiee! Gemax 
levels of factors. The values of these parameters are presented in Table 1. 


Table 1 
Values of experiment factor levels 
L level U level (+1 
Factors Code a Main level (0) pper level (+1) 
f,, mm/tooth x 0.04 0.08 0.12 
V., m/min X2 80 90 100 
ae, mm x3 1 2 3 
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The experiment was planned according to the non-composite design proposed by Box and Behnken. This plan is a 
sampling of rows from a complete 3‘ type factor experiment. It includes 15 experiments. The matrix of the plan is 
shown in Table 2. 


Table 2 
Experiment Design Results Matrix 

No Factors Output parameter 
Xo xy Xo x3 Hirxe | res | eerws | wa? | we | ws” Ra, pm 
1 | 1 -1 -1 0 1 0 0 1 1/0 0.9487 
2 | 1 1 -1 0 -1 0 0 1 1 | 0 1.5583 
3} 1 -1 1 0 -1 0 0 1 1 | 0 1.2183 
4 }1 1 1 0 1 0 0 1 1 | 0 1.7555 
5 | 1 -1 0 -1 0 1 0 1/0; 1 1.5674 
6 | 1 1 0 -1 0 -1 0 1; 0; 1 1.9084 
7 | 1 -1 0 1 0 -1 0 1/0; 1 1.5198 
8 | 1 1 0 1 0 1 0 1);0/; 1 1.9505 
9 | 1 0 -1 -1 0 0 1 0} 1 1 1.8912 
10 | 1 0 1 -1 0 0 -1 0 |} 1 1 1.0914 
11) 1 0 -1 1 0 0 -1 0} 1 1 1.6132 
12 | 1 0 1 1 0 0 1 0} 1 1 1.7664 
13 | 1 0 0 0 0 0 0 0}; 0] 0 1.5668 
14) 1 0 0 0 0 0 0 0/0; 0 1.5081 
15 | 1 0 0 0 0 0 0 0/;0;] 0 1.5046 


The coefficient vector is determined from the formula: 
B= (XTX) UX'TY. 
Variances S? (b;) of the regression coefficient are found from the expression: 
S?(bj) = cS}, 
where c;; — diagonal elements of matrix XX"'; S’y — empirical regression variance. 
The average value and variance are determined by the vector of experiments at the center of the plan. 
Average value: 


—_ Lheo Yas; 

= 

where Y,, — vector of response function values in the center of the plan. 
Dispersion: 


Lieou,;-y)? 
2 _ 4i=0 Lay. 
3 = 


Based on the results of experiments carried out according to the plan (Table 2), it is possible to determine the 


constant coefficients of the equation: 
Y= Go Xo + 91° X1 + Go X2 + G3°Xz + Giz" %1*X2 + Gig X1*X3 + Gog" X2°X3 


+911 °° X2 + G22 XZ + 933° X53 


no 
_i >» 
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u=1 
N 
GHA > XipVj> 
jaa 


N 
Ji = D >: XijX1jVj > 


j=1 


Machine building and machine science 


103 


http://vestnik-donstu.ru 


Advanced Engineering Research 2022. V. 22, no. 2. P. 99-106. ISSN 2687-1653 


no 


N k N 
2 2 1 
93 =B > XY) +C > > XiGVji — > You 
= Pp* No 
u=1 U 


i=1 j=1 =1 

where ny — number of experiments in the center of the plan; uw — number of parallel experiments in the center of 
the plan; yo,— value of the response function in the u-th experiment; N — number of experiments in the planning 
matrix; j — number of experiments in the planning matrix; i, /— factor numbers; x;, x,; — coded values of the i-th and 
I-the number of factors in the j-th experiment; y; — value of the response function in the j-th experiment. 

It follows from this: gg = 1.388; g, = 0.24; gz = —0.022; 93 = 0.268; gy, = 2.413; go2 = 2.267; g33 = 3.072. 

Thus, the model of the dependence of the deviation of the surface roughness on the processing parameters and the 
cutting width will have the form: 

Ra = 1.388 + 2.413 - f? — 0.24: f. 


Research Results. The dependences of the surface roughness on the cutting modes are shown in Figure 3. 


m>2 
Mm <2 
GB <is 
G@<i6 
@<i14 
Mm <12 


a) b) c) 


Fig. 3. Dependences of the deviation of roughness parameter R,: a — on feed to tooth f, and on cutting speed V,; 


b — on cutting speed V, and milling width a,; c — on milling width a, and feed to tooth f, 


The analysis of the obtained results shows that the dependence of roughness on feed has a linear character over the 
entire investigated range of cutting modes (Fig. 3 a, b). At the same time, with an increase in feed, the roughness 
increases proportionally, which is consistent with the basic concepts of the cutting theory. The maximum value of the 
roughness parameter R, is 1.85 um, which corresponds to the most common requirements for finishing machining in 
mechanical engineering. 

The dependence of roughness on the cutting speed V, (Fig.3.a, b) has a parabolic character with maximum 
roughness values in the range of cutting speeds of 90-95 m/min. This can be explained by the fact that in the specified 
range of cutting speeds, there is an effect of built-up material on the cutting edge. To verify this hypothesis, it is 
required to conduct additional studies with temperature control in the cutting zone and the use of a high-speed video 
camera to monitor the process of chip formation and possible build-up formation. 

The impact of milling width a, on the surface roughness is also parabolic in nature (Fig. 3 b, c), but the maximum 
values of roughness are reached at a, = 0.25 mm and a, = 3 mm, and the minimum value is at a, = 1.8—2.2 mm. It can 
be assumed that in the specified cutting width range, the treatment process is less affected by regenerative vibrations, 
which also depend on the design and geometric parameters of the cutting tool. 

As a result, it can be stated that with the established cutting modes, it was possible to achieve a surface roughness no 
higher than R, = 1.85 without loss of process performance. 

Discussion and Conclusions. A mathematical model has been developed reflecting the effect of cutting modes on 


the surface roughness through turn-milling of high-speed steel with a monolithic cutter. The results of the conducted 
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research can be used to determine the optimal cutting conditions that provide a given surface quality in the manufacture 
of real parts under production conditions. 

It is recommended to conduct additional studies with temperature control in the cutting area and vibrations. The 
assignment of cutting modes should be carried out on the basis of the results of the modal analysis of the technological 
system to reduce the effect of regenerative self-oscillations on the roughness of the treated surface. 
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Abstract 

Introduction. Over the past decade, global demand for pressure vessels has increased significantly, specifically in such 
industries as aviation, space, chemical, and oil and gas. Being under the constant impact of high internal pressure, the 
walls of the tanks are under increased stress, which can cause their sudden destruction. To eliminate this probability and 
improve the strength characteristics, the tanks are made in the form of metal cylinders with an internal coating of 
composite material consisting of resin reinforced with carbon fibers. This article aimed at studying the effect of the 
angle of inclination of carbon fiber on cylindrical tanks and determining the maximum destructive pressure using the 
finite element method of ANSYS program. 

Materials and Methods. Using the ANSYS program, a finite element model of a tank was created. It has a central part, 
which is a metal cylinder with an internal coating of composite material consisting of polymer reinforced with carbon 
fibers. At the ends of the tank, spiral wound hemispheres were placed. In these studies, SHELL 99 was used to model 
the layered composite material. The Tsai-Wu theory was used to determine the pressure tank failure criterion. 

Results. The cylindrical tank model was calculated for two types of fiber winding paths: annular and spiral, at different 
angles of their inclination. The results of the pressure value analysis for different fiber inclination angles showed that, 
starting from the angle value of 0° and up to 45°, it increased, and then, up to the angle value of 65°, it began to 
decrease. The critical pressure value for a carbon fiber reinforced tank was 207 MPa, which was obtained at a fiber 
angle of 45°. 

Discussion and Conclusion. Analysis of the studies showed that at a fiber inclination angle of 45°, the value of the 
maximum stress turned out to be the smallest, and the maximum possible destructive pressure at the same angle was 
207 MPa. It follows, that the optimal fiber orientation angle to provide safe operation of the high-pressure tank is + 45°, 
and the carbon fiber tank, calculated at the same fiber winding angle, has the maximum strength value. 

Keywords: high-pressure tank, computer model, winding angle, composite coating, carbon fiber, polymer binder. 
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Introduction. The use of high-pressure carbon fiber vessels has found widespread use in various industries due to 
its distinctive properties, such as low weight and high strength. Therefore, recently the demand for such vessels has 
increased significantly in cases where weight plays a major role [1, 2]. 

The critical areas of application of this type of tanks are aerospace, aviation, and chemical engineering. In addition, 
fiber-reinforced tanks are widely used for increase in pressure during the transportation of oil and gas. In many cases of 
such use, tanks are exposed to high internal pressure, which may result in a step-like increase in stress on the vessel 
walls and their sudden collapse, causing great damage to material and human resources [3, 4]. 

In [5], a number of studies were carried out aimed only at exploring a pressure tank made of multilayer composite, 
during which the expected collapse resistance was determined. In [6], the behavior of a rotating composite pressure 
vessel under the internal pressure and axial load was studied. In [7], the effect of thermal loads on a multilayer 
composite pressure tank was described. In [8], the behavior of a polygonal composite pressure tank of five different 
shapes under the action of internal pressure of various modes was studied; and in [9], the performance of a composite 
pressure tank was explored under the influence of transverse loads. 

The design of a composite tank is a challenge; therefore, priority factors should be selected to conduct a complete 
and accurate analysis. We studied a pressure vessel reinforced with several layers of carbon fibers and subjected to 
internal pressure loading. To determine the maximum stresses and displacements at operating pressure, to identify the 
expected limiting pressure causing destruction, as well as to identify the optimal angle of orientation of the fibers, a 
finite element model of the tank was created using ANSYS customized application. 

Materials and Methods. Theoretical Research. The high-pressure tank made of carbon fibre reinforced 
plastic (CFRP) consists of a central part, which is a metal cylinder with an internal coating made from a polymer 


reinforced with carbon fibers. The end surfaces of the tank have the shape of hemispheres with spiral winding (Fig. 1). 
ee 
See 
Fig. 1. Longitudinal section of the tank 

Tanks operating at high pressures and reinforced with carbon fiber are manufactured by the method of filament 
winding (Fig. 2). To obtain the required reinforcement stability, the fibers are sent to a moving trolley with careful 
selection of their coordination, and then wound onto a cylindrical surface. The stability of fiber coordination is 
influenced by several factors: temperature, surface shape and treatment, as well as the degree of adhesion of fibers to 
the matrix. The winding angle is controlled by the speed of the trolley and the rotation speed of the cylindrical drum. To 


obtain high operational and strength properties, the inner cylindrical surface of the tank is covered with several layers of 


fibers [10]. 


Imad Rizakalla Antipas, et al.Using the Finite Element Method to Simulate a Carbon Fiber Reinforced Polymer Pressure Vessel 
a 


Fiber angle 


Fig. 2. Fiber winding scheme 


The angle of the fibers has a significant effect on the properties of the vessels; therefore, to find the appropriate 
angle for each part of the vessel is of critical importance. The fiber orientation angle is determined by the required 
amount of friction between the fibers and the composite material layer, as shown in the relation: 


a(R R-R, ) 
a(R)=sin (4) +0( i ) ; (1) 


where R — distance between the center and the point of the layer; R, — central axis radius; R,, 


— radius on the 


tangent to the surface of the cylinder at 6 =0 [8]. 

Tsai-Wu Destruction Criteria. When studying and simulating a high-pressure tank made of composite material, 
the destruction criteria according to the Tsai-Wu theory [11] were used. The implementation of equation (1) is required 
to study the expected fracture of orthotropic materials according to the Tsai-Wu theory: 

Fo, + Foy + Feo, FiO ee +2F,,0);0y =1. 


Elastic properties are determined by four independent constants: E,,, Ex, Gi2, Vi2, presented in Table 1. 


Table | 
Characteristics of LY5052/T300 composite material 
Properties Carbon fiber T300 Epoxide LY5052 
Elasticity modulus 230 GPa 3.0 GPa 
Tensile Strength 3.5 GPa 71.0 GPa 
Density 1,760 g/cm" 1.14 g/cm? 
Forces are calculated from the following equations: 
1 1 1 1 1 
ss -aS Se 
x, x, Y, Y, X, x, 
1 1 1 
Fy 7 FF» » 
where X, — tensile force in the longitudinal direction; Y, — tensile force in the transverse direction; X, — pressure 


force in the longitudinal direction; Y, — pressure force in the transverse direction; § — shear force. 
The maximum destructive stress is reached when one of the following ratios is met [12]: 


O; 0. 
1>1; 2 


cy 


>1: Ti. >I. 
S 


Material properties and finite element modeling. The composite material used for the tank under study is 
polymer reinforced with T300-type carbon fiber, and LY5052 epoxy is used as polymer material. Composite materials 
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are orthotropic in nature; therefore, the process of modeling them with finite elements is more complicated than 
isotropic materials, such as aluminum and steel. 

Figure 3 shows a finite element model of a high-pressure tank made of carbon fiber, whose inner layer consists of an 
aluminum alloy reinforced with eight layers of Carbon/Epoxy T300/LY5052 composite material. 

The model has the following dimensions: 

— tank length —1,200 mm; 

— tank diameter in the center — 300 mm; 

— total thickness — 64 mm; 

— thickness of one layer — 6.5 mm; 

— lining thickness — 0.12 mm. 


Fig. 3. Finite element model of the CFRP tank 


To study high-pressure tanks made of composite materials, it is very important to select the appropriate type of finite 
elements. ANSYS program contains SHELL and SOLID finite elements required for modeling layered composite 
materials. For this study, SHELL 99 program was used, which accelerated the calculation of a structure with up to 
250 layers. This element is a multi-level linear structure with eight nodes and six degrees of freedom. It allows the user 
to determine the flexibility, slope of layers, and density of each layer. 

When studying composite materials, the formation of layers is one of the major issues, since each layer has its own 
angle of inclination, and the fibers in each layer have different angles of inclination; therefore, the properties of each 
layer should be determined separately. 

The formation of layers requires the characteristics of the material, the number of layers, the fiber inclination angles, 
the thickness of the layer, and the number of integration points in each layer. Here is some comparative information 
about the properties of the materials used: CFRP and aluminum 6061. Density: CFRP — 1,570 kg/m’, aluminum — 
2,750 kg/m’. Table 2 shows mechanical characteristics of the materials used: elastic modulus E; shear strength G; 


tensile strength V. 
Table 2 


Mechanical properties of CFRP and aluminum 6061 


Parameters in eae of coordinate CFRP eminanrebel 

Ex 128 7070 

Ey 10.5 70 

Gyy 5 70 

Gyz 5 

Gzy 5 

Vyy 0.27 0.3 

Vyz 0.4 0.3 

Vex 0.02 0.3 


[ +25° /-25° |, [ +30° /-30° |, [ +35° /-35° |, [ +40° /-40° |, [ +45° /-45° |, [ +50° /-50° |, [ +55° /-55° |. 
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Figure 4 shows a sequence of layers with the fiber inclination at an angle of + 45°, implemented in ANSYS 


program. 


AN 


Jul3 2014 


LAYER STACKING 


ELEM ° 0 21:24:22 
SECT - 1 
LAYERS - 
TOTAL .- O 
SHOMM 1 TO 
FROM 0 
Layers Materials 
7 7 7 7 “? 1 
1 , 
There a 
45 
-45 


45 


Fig. 4. Sequence of layers with the fiber inclination at an angle of + 45° 


Digital simulation. Calculation of stresses and displacements. The operation of the high-pressure tank was 
analyzed using the destruction criteria according to Tsai-Wu theory. Internal working pressure of 35 MPa was used to 
calculate the maximum stresses and displacements. The greatest stress of the tank was observed at the angle of 
inclination of the annular fibers — 0° and various angles of inclination of the spiral fibers: 


[ +25° /-25° |, [ +30° /-30° |, [ +35° /-35° |, [ +40" /-40° |, [ +45° 45°, [ +50° /-50° |, [ +55° /-55° |. 


Figure 5 shows the distribution of equivalent stresses and displacements in the tank at the fiber inclination of + 45°. 


Tables 3, 4 show the maximum and minimum values of displacements and equivalent stresses in the directions of 
X, Y, Z coordinate axes. 


‘ a rc 
MODAL SOLUTION AN | MODAL SOLUTION AN S 
STEP=1 Jul 3 2014 STEP=1 Jul 3.2014 
SUB=1 13:01:07 SUB=1 13:00:15 
TIME=1 TIME=1 
USUM SEQV (AVG) 
(AVG) DMX=6.473 
RSYS=0 SMN=11.886 
DMX=6.473 SMX=463.756 
SMX=6.473 
0 1.439 2.877 4,316 5.754 11,886 112.302 212.717 313.133 413.549 
719265 2.158 3.596 5.035 6.473 62.094 162.51 262.925 363,341 463.756 
a) b) 


Fig. 5. Equivalent distribution in the tank: a — stresses; b — displacements 
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Table 3 
Maximum and minimum displacement values 
in the directions of X, Y, Z axes 
Direction of deformation Minimum, mm Maximum, mm 
total vector 0 6.473 
along X-axis -3.914 3.914 
along Y-axis -6.473 0 
along Z-axis -3.906 3.909 
Table 4 
Maximum and minimum stress values in the direction of X, Y, Z axes 
Direction of deformation Minimum, MPa Maximum, MPa 
total vector 11.886 463.756 
along X-axis -81.852 528.485 
along Y-axis -155.587 341.607 
along Z-axis -80.307 530.07 


Determination of the expected destructive pressure and the optimal fiber inclination angle. The loading of the 
tank was carried out through gradually increasing the internal pressure, starting from its operating value — 35 MPa. 
Then, the obtained design maximum internal pressure stress was compared to the final allowable stress under the 
condition: o,,,, <o, , where: o,,,, — design maximum stress; o,, — final allowable stress. The permissible stress for 
CFRP pressure tanks is 1,210 MPa. 

The cylindrical tank model was calculated for two types of fiber winding paths: annular and spiral, as well as at 
different fiber inclination angles: 


[ +25° /-25° ], [+30° /-30° ], [+35° /-35° ], [+40°/-40° ], [+45° /-45° ], [+50° /-50° ], [ +55° /-55° ]. 
With an annular winding path, the fiber inclination angle to the axis of the cylinder was 0 °. Figure 6 shows a finite 
element model of a cylindrical high-pressure tank. 


ELEMENT 1 


Fig. 6. Finite-element model of a cylindrical high-pressure tank 


Research Results. Figure 7 shows the change in the maximum stress; starting from the angle value from 0° to 45°, 
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it decreased, and then increased to the value of 65°. It was found that the maximum stress in a cylindrical CFRP 


pressure tank was the lowest at the fiber inclination angle of + 45°. 
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Fig. 7. Dependence of maximum stress on the fiber inclination angle 


Figure 8 shows the results of the analysis of the pressure value for various fiber inclination angles: starting from the 
value of angle 0° and up to 45°, the pressure increased, and then it decreased to the value of angle 65°. The maximum 
pressure that a CFRP tank can withstand without destruction is 207 MPa. This pressure occurs at the fiber inclination 
angle of + 45°, i.e., it is the expected value of the destructive pressure for the tank. 
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Fig. 8. Dependence of pressure on the fiber inclination angle 


Figure 9 shows the distribution of stresses in the tank with fiber inclination angle of + 45°, and the maximum 
pressure stress is 1,213 MPa, which is higher than the allowable stress for CFRP pressure tanks 1,210 MPa. 


—— 
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Fig. 9. Stress distribution in the tank under pressurization and at fiber angle of + 45° 113 
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Discussion and Conclusions. A finite element model of a carbon fiber reinforced tank was designed and analyzed 
using ANSYS program with the application of SHELL 99 element in the course of the simulation process. Several 
models were created with different fiber angles. The maximum stress and collapse pressure for the high-pressure vessel 
were calculated using the Tsai-Wu criteria. It was found that the maximum stress was the smallest at the fiber 
inclination angle of + 45°, and the maximum possible collapse pressure under the same conditions was 207 MPa. This 
indicates that the optimal fiber angle for safe operation of a pressure vessel is + 45°, and a CFRP vessel rated at the 


same fiber winding angle has maximum strength. 
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Abstract 
Introduction. Domestic and foreign works on the milling of complex-profile surfaces with a ball-end tool were 


analyzed. Methods of surface quality control and ways to provide amplitude parameters of roughness, based on research 
data and field experiments, were considered. Theoretical provisions on the determination of cutting forces and the 
results of vibroacoustic diagnostics were presented. 

Materials and Methods. The methods of correlation analysis, comparison and generalization of the results were 
applied. The data were calculated at different tool angles, taking into account the instantaneous cutting forces, and were 
fixed in the range of values of the variable feed per tooth (fz) and the angle of inclination of the surface (y). The 
vibroacoustic diagnostic data and theoretical data of the presented model at different tool inclination angles were 
verified by experiment. Consequently, such methods can be used to predict surface roughness parameters. 

Results. The relationship between cutting forces, tool inclination angle, and vibroacoustic diagnostics data was found. A 
model of the cutting force and tool displacements was formulated taking into account the inclination of the surface. The 
optimal range of the inclination angle of the tool to the surface to be machined, at which the minimum values of the 
amplitude parameters of roughness were achieved, was determined. The sound vibrations obtained empirically, 
presented in spectral and wave forms, were in good agreement with data from other sources. This allowed us to 
conclude about the feasibility of forecasting and monitoring roughness parameters in real time through acoustics. 
Discussion and Conclusions. It was established that the growth of forces in the direction a,(X) and fz(Y) was observed 
at y > 40°. This was due to the distribution of the components of the cutting force along the cutting edge and depended 
on the inclination of the surface. The amplitude parameters decreased when the angle increased from 10 to 40 degrees. 
The found interrelations of force analysis, processing directions, and vibroacoustic diagnostics have validated the use of 
vibroacoustic diagnostics to predict surface roughness. Acoustic diagnostics, regardless of the layout of technological 


equipment, enables to quickly adjust the sound device and assess the impact of cutting modes on roughness parameters. 
Keywords: surface roughness, ball-end tool, vibroacoustic diagnostics, milling, tool inclination, vibration, cutting forces. 
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Introduction. Specialists do not have enough software tools that could help predicting the value of the cutting force, 
vibration amplitude, etc., when programming machines with numerical control (CNC). As a result, it is required to 
select the feed rate based on experience, setting a constant feed rate for the entire processed complex surface. This can 
reduce the processing efficiency and surface quality. To solve this problem, it is required to build a highly accurate and 
reliable model for predicting surface roughness. 

Currently, two main approaches are used in such forecasting. In the first case, photoelectric sensors transmit 
roughness data online. This method efficiency is quite low because chips are formed during the cutting process, and a 
lubricating fluid (lubricoolant) is supplied to the area of the recorded microrelief. The second group includes control 
methods using artificial intelligence, as well as theoretical modeling and modeling of empirical regression analysis. In 
practice, methods from the second group are most common. 

The authors of paper [1] explored the spherical tool milling of aerospace components made of aluminum alloy LM6 
and developed a model using ARMAX structures (autoregressive moving average with exogenous inputs) to predict 
surface roughness under various processing modes. In [2], a mathematical model of the trajectory of a spherocylindrical 
tool for processing mushroom-shaped blades of a steam turbine was presented to improve the processing efficiency and 
surface quality. The authors [3] created an analytical model of the milling with a spherocylindrical tool and considered 
the situation of vibration occurrence to predict the limits of tool stability at different angles of inclination and spindle 
rotation speed. 

In [4], the processing of the impeller on a five-axis CNC machine using a spherical cylindrical tool was studied. The 
tool trajectory and optimal cutting modes had to be developed due to the complex curved surface and the close location 
of the wheel blades, since in this case, it was impossible to use standard CAM systems (computer-aided manufacturing - 
an automated system or its module for the preparation of control production programs). 

In [5], a system for optimizing the feed rate for three-axis milling with a spherical tool was described. The authors 
used the geometric modeling capabilities of B-rep ACIS: miter cutting was converted to an orthogonal model for 
various tools and workpieces. The authors [6] created a mathematical model of a constant cutting force [7] for 
processing a complex curved surface with a spherocylindrical tool on a CNC machine through maintaining a constant 
cutting speed at different tool inclination angles. Paper [8] presents a dynamic model of the milling with a spherical 
cylindrical tool, including cutting modes and tool edge wear. The authors [9] developed a real-time monitoring system 
with error compensation to improve accuracy in the production of free-form parts. In [10], the energy consumption of 
various cutting strategies was studied to evaluate the efficiency of each cutting strategy. The authors of this work have 
established that the optimal strategy enables to provide the specified parameters of surface quality. Many scientists used 
the Taguchi method to optimize processing parameters [11—13]. Some of the surveys are narrowly focused, and more 
work is needed to scale up the results [14-18]. 

The analysis of literature sources found the absence of recommendations on using the method of forecasting and 
monitoring roughness in real time. This is even more so in cases where sound vibrations in the air are taken into 
account. 


Materials and Methods. To determine the cutting forces, a model from [7, 19] is used: 


dF =Kydl, + KAS, 
dF, =Kyedl; +KyedSz, 


(1) 
dF, =Kgedl; + KgcAS.;. 


Here, Kj, Kye, Kae — edge specific coefficients; K,,, K;-, Kae — shear specific coefficients; dl; — incremental cutting- 
edge length /; [20]; S.; — cross-sectional area of the cut segment. 

In this model, a set of curvilinear coordinate systems, perpendicular to the tangents to the spherical surface, is used 
to set the resulting force acting on the i-th infinitesimal segment of the cutting edge. 


The geometry of the spherocylindrical end mill is schematically shown in Figure 1. 
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Fig. 1. Diagram for determining: a — cutting forces of a spherocylindrical cutter; 


b — segment to be cut; c — cross-sectional area of the given figure [16] 


The diagram shows the angle of inclination 4 of a milling cutter with a spherocylindrical initial surface of the tool 
along the cutting edge. Here, we also see: 

— tangential dF, radial dF, axial dF; cutting forces acting on the j-th tooth, and the coordinates of the tool for the 
spherocylindrical end mill (Fig. 1 a); 

— chord L,,;,, segment angle ®,, arc length /,,, (Fig. 1 b); 

— cross-section area (Fig. | b, c). 
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= . ht 2 
Sj = Lyon * Maver: @) 


where L>,;, — length of the second chord, mm; h’,,,., — average chip thickness, mm. 


(3) 
I very = «SZ + SIMA. (4) 
Length of the first chord L;,;, arc length /,,,, and circle segment angle ®;: 
L, = 2-sin- lwe . R, (5) 
ch 2R 
®O,-R 
lore = > 
ze (6) 
@, = 2arccos [1 -2- =). 
2R (7) 
Angle A (Fig. | c) for calculating the average chip thickness: 
A = arctan [4] : 
Lin / 2 (8) 
Radius vector r of the current point of the initial tool surface: 
Rsing cosQ. 
Rsing sinQ 
ra aly (9) 
Reos@ 
1 


The equation of the cutting edge is solved as the equation of a curve located on the original instrumental surface and 
crossing the meridians at angle 190°. A parametric equation of the form © = () describes a certain line on the 
sphere. Cosines of the angles between the axes of XYZ coordinate system of the spherocylindrical cutter and the tangent 
line to the curve Q = Q (q) on the initial tool surface: 

dX cos@ cosQ do — sing sinQ do 


cose =, 
dS Jde + sin’ dQ? 


dY cos@ sinQ d@ — sing cosQ do 


cosB = a P 
dS Vdo + sin dQ? (10) 


dZ — sing d@ 


cos6 = = 7 
dS do + sin“@p dQ 


where dS — cutting edge arc differential. 
Cosines of the angles between the coordinate axes and the tangent to the meridian (at Q = const) on the initial tool 


surface are represented as Radzevich’s dependences [21]: 


dX 
cosa = —— = cos@ cosQ, 
dS 


cospB = — = cos@ sinQ, (11) 


dZ ; 
cos6 = — = -— sing. 


ds 
We set the condition: Q = i. Then, the cutting-edge equations in parametric form: 
Reosi 
chq (2% + C)- 
Rsinx (12) 
chq (A + C) 


Z = Rthg(r% + C), 
where ch — hyperbolic cosine, th — hyperbolic tangent. 
To determine the infinitesimal length of the cutting edge, we can use the expression proposed in [7, 16, 19]: 
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dl = sr’ (w) + (r' (wy) + (<' (wy) ay, (13) 
where 
r(y) = Vl - (y cot a — 1) 
toe —R (y cot A — 1) cotr . (14) 
V1 - (wy cot A — 1) 
z’(w) = X* + ¥7 + (R- ZY cota 
Main Text 


Experiment Details. At the preparatory stage, the literature was studied [22-27]. The data for planning the 
experiment are described below. Machining was carried out in passing and counter directions (Fig. 2), on blanks with 
the properties of steel 45 (foreign analogues — C45, 1045). We used a spherocylindrical end carbide milling cutter with 
a diameter of D = 8 mm, with two teeth from Sandvik Coromant (z = 2). The ratio of the overhang of the tool installed 
in the chuck was taken as | / D, equal to 4. The milling tool made of fine-grained tungsten carbide had anti-wear coating 
TiAIN and parameter 4 = 30°. Feed per tooth fz = 0.2 mm/tooth, allowance for all samples a, = H = 0.2 mm, lateral 
pitch a, = 0.2 mm. 

To provide equal cutting speed V, at different angles of inclination (y) of the tool to the surface to be processed, the 
rotation frequency (n) varied in the range from 1478 min ' to 8000 min‘. 


Z 
a, 
~ el 
a) 


Fig. 2. Types of the conditional contact surface of workpiece and tool at different values of tool inclination angles in the direction: 


b) 


a — passing; b — counter 


Angle y was selected from the range from 10° to 50°; therefore, the calculation of the effective diameter Deap (Lich) 
was carried out using the formulas: 


= 2Dcosm = 2. R - cos (90° — y) 
y > O° 40r (15) 
l 
D 7 = 2-sin:-—.-R. 
cap Lich OR 

In the experiments, a DMG DMU 50 CNC machine (Heidenhain control system) was used. The surface roughness 
after treatment was measured using profilometer-profilograph-contorograph Surfcam 1800D. In accordance with the 
current standards', measurements were carried out along various routes. To find the maximum values of the surface 
roughness parameters, 50 % Gaussian filtration was used. Vibroacoustic diagnostics was performed according to GOST 
R ISO 7919-3-997 using spectrum analyzer Zet 017-U2. As an output evaluation of the processing efficiency, amplitude 

A of the vibrations of the spindle assembly (mm), varying with time (tf, s) was used. 


' GOST 27964-88. Measurement of surface roughness parameters. Terms and definitions. Moscow, 2018. 13 p. (In Russ.) 
? GOST R ISO 7919-3-99. Mechanical vibration of non-reciprocating machines. Measurements on rotating shafts and evaluation criteria. Moscow, 


2000. 12 p. (In Russ.) 
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Research Results. The experiment has shown that the real surface roughness differs significantly from the data of 
theoretical models — the results of the kinematic-geometric projection of the tool on the workpiece. Currently, there is 
a growing demand for details of a spatially complex configuration and surface quality requirements. In this regard, it is 
worth focusing on the requirements of predictability of machining results, expressed, among other things, by roughness 
indicators. 

Tables 1 and 2 present the results of the experiment of amplitude parameters of roughness described in 
GOST R ISO 4288”. 


Table 1 
Roughness parameters depending on the machining direction 
with a spherocylindrical tool under climb milling 
Amplitude parameters of roughness, um 
Angle, ° Ra Rq Rz Rt Rp Rv 
10 0.457 0.569 2.912 5.000 1.671 1.240 
20 0.597 0.735 3.936 5.159 2.335 1.600 
30 0.449 0.561 3.007 4.960 1.808 1.200 
40 0.479 0.590 3.232 4.239 1.823 1.407 
50 0.787 1.027 6.615 12.68 4.119 2.496 
Table 2 
Roughness parameters depending on the machining direction 
with a spherocylindrical tool under up-cut milling 
Amplitude parameters of roughness, 1m 
Angle, ° Ra Rq Rz Rt Rp Rv 
10 0.806 1.006 4.776 5.960 2.776 2.000 
20 0.752 0.935 5.207 6.440 2.847 2.359 
30 0.509 0.647 3.727 5.239 2.223 1.503 
40 0.628 0.805 4.832 7.400 2.912 1.919 
50 0.730 0.915 5.111 7.639 3.039 2.072 


Analyzing the data of Tables | and 2, we note that the prevailing influence on the “height” of the roughness is 
exerted by the cutting direction. The quantitative value of the amplitude parameters is “lower” under climb milling. 
Hereafter, we will talk only about this type of milling. Based on the data in Table 1, a correlation analysis was 
performed, and a correlogram was made up (Fig. 3). The found dependences were correlated with the results of 
papers [28, 29]. 


— 0.9...1.00 0.8...0.89 --~  0.7...0.79 


Fig. 3. Correlogram of the amplitude roughness parameters after milling under climb milling 


The equations obtained during the correlation analysis helped us to move away from the direct normalization of 
roughness parameters. At the same time, analytical dependences provided minimizing the number of controlled factors. 
This reduced the complexity of the transition to real-time monitoring of sound vibrations. 


3GOST R ISO 4287-2014. Geometrical Product Specifications (GPS). Surface texture. Profile method. Terms, definitions and surface texture 
parameters. Moscow, 2015. 20 p. (In Russ.) 
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In the presented study, the cutting forces (Fx, Fy, Fz) were calculated (Fig. 4). It can be seen that the inclination 
angle of the surface y affected significantly the cutting forces. 


Fi,N Fi, N 
100 100 
—Fy Fy 
50 Fz 50 —F, 
Fx Fx 
0 0 
—50 50 
—~100 i 100 ts 
0 0.005 0.01 0.015 0.02 t, 0 0.005 0.01 0.015 0.02 
a) b) 


Fig. 4. Cutting forces calculated for: a— y = 10; b— y = 40 


The data obtained allowed us to establish that the forces decreased with an increase in the tool inclination angle 
from 10° to 40°. In addition, force Fx was most sensitive when the inclination angle of the surface changed. It should 
also be noted that the cutting forces calculated according to the above-mentioned method were in good agreement with 
the empirical results obtained [6, 17, 24, 26, 27, 30]. 

With an increase in angle y, the maximum area of plunging of the tool into the workpiece decreased (determined by 
the working angle and the active number of teeth z,). Consequently, when finishing milling with a spherocylindrical 
tool, “pulsating” forces may be present, since for the inclination angle of the surface y > 0°, the number of active teeth 
is often less than one (z, < 1). 

To establish z,, sound vibrations during cutting were analyzed. The authors of [31] proposed to divide the spectrum 
of vibrations of working bodies into ranges: 

— 20...300 Hz — low frequency, 

— 300...1500 Hz — medium frequency, 

— 1500 Hz and more — high frequency. 

The experimentally obtained range of sound vibrations is shown in Figure 5. 


Fig. 5. Experimentally obtained spectral type of sound wave (segment is | sec) 


The vibrations of the mounting elements and feed drives should be attributed to the low frequency range, the 
frequency of tool, spindle, etc. — to the medium and high frequency ranges. In this work, from the point of view of 
acoustic diagnostics, optimal cutting modes were determined that supported the stable operation of the technological 


system. 
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The frequency corresponding to the cutting frequency was pre-determined (Fig. 6). Filtering was applied to sound 
vibrations under: 

— spindle assembly movement (Y, Z), 

— rotation (n), 

— table movement (X), 


— noise in the laboratory. 


Fig. 6. Spectrogram of the sound wave after filtering (segment is equal to 1 sec) 


Some spread in values (Fig. 6), associated with a huge number of factors that were difficult to exclude when 
conducting a full-scale experiment, is noteworthy. 

The middle of the pronounced frequency range shown in Figure 5 corresponded to a rotation of 1500 rpm with a tool 
tilt of 40°. Acoustic diagnostics confirmed that the cutting went on almost continuously, i.e., without shock “pulsating” 


loads. The latter situation is shown in Figure 7. Here, the allocated evaluation time is 0.5 s, the number of periods is 25. 


Fig. 7. Filtered fragment of the sound wave based on the sound of cutting at y = 40° (segment is 1/2 sec) 


Figure 7 shows an ambiguously interpreted sound wave. However, it is quite easy to characterize it if we delve into 
the essence of mechanical milling. Under the climb milling and rounding of the cutting edge, plastic deformations of the 


metal with residual stresses remain on the treated surface. Even with proper coolant supply, temperature deformations 
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occur in the surface layers of the material and distort the crystal lattice (residual stresses of the HI kind). Along with 
this, after removing the tool, the stretched upper layers of the metal acquire residual compression stresses. If we 
describe the machining from the kinematics of the process, of note, the chip removal factor, namely value ©, — the 


cutting angle of the material with the inclined spindle in the passing direction (Fig. 8). 123 


http://vestnik-donstu.ru 


Advanced Engineering Research 2022. V. 22, no. 2. P. 116-129. ISSN 2687-1653 


Fig. 8. Segment 2 to be cut off during machining: a — top view; b — projection of the segment being cut [16] 


We specify that the tool cut in zone 2. The data obtained were validated by the simulation results, but were 
considered for the ideal case. 

Along with these factors, theoretical and empirical discrepancies were often caused by vibrations associated with 
tool deflection and wear, as well as with the geometry of the cutting zone. Figure 7 shows the process of cutting the 
workpiece material with the j-tooth of the tool, which, along with other factors, determines such an interesting shape of 
the resulting wave. 

In addition to the above, Figure 9 presents dependences which show that with a decrease in the amplitude of 


vibrations, the angle of inclination increased from 10° to 40°. 
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Fig. 9. Vibration amplitudes depending on the surface inclination angle: solid line — 10°, dotted line — 40° 


M. R. Gimadeev and V. V. Gusliakov* analyzed the international practice of profile assessment of surface 
roughness. The authors of work [32] put forward proposals for the transition to a three-dimensional assessment. In the 
presented study, roughness measurements were carried out in an orthogonal coordinate system to make a 
comprehensive assessment. 

The analysis of the profilograms obtained (Fig. 10) allowed us to conclude about the regularity of the surface 
microrelief after processing with a spherocylindrical end mill. In this case, it is impractical to take into account the 
roughness of the tool, since its manifestations are insignificant. Some additional increase in roughness was observed in 
the direction of the side pitch a,, which verified the increase in elastic squeezes of the tool. The presented profiles, 


where the profile periods (peak-to-valley) correspond to a given feed value fz = 0.2 mm/tooth, testify to this. 


4 Gimadeev MR. Analiz zarubezhnogo opyta pri frezerovanii sferotsilindricheskim instrumentom. Khabarovsk: Informatsionnye tekhnologii 
XXI veka; 2019. P. 330-334. (In Russ.) 
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Fig. 10. Profilograms of surfaces after milling at different inclination angles y: 
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a— y= 10; direction f,; b — y = 10; direction a,; c — y = 40; direction f,; d—— y = 40; direction a, 


Photomicrographic images of the treated surface have an obvious hexagonal shape (Fig. 11). The quality of the 
machined surface in the direction of the side pitch a, may decrease due to the cutting tool retraction and an increase in 
the amplitude of the technological system vibrations under milling. This is validated by the results of vibroacoustic 
diagnostics. 125 
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b) 
Fig. 11. Milling scheme and photo of the treated surface at: a— y = 10; b— y =40 


We emphasize once again that the experimental results correlate with the data obtained by many authors, and are 
consistent with the theories of elasticity, shaping, material cutting, and mechanical engineering technology. 

Discussion and Conclusions. The influence of the inclination angle of the treated surface on the cutting forces, 
vibrations and microrelief parameters was studied. A model of cutting force and vibration was created, including 
kinematic and geometric parameters. The study has shown that quantitatively and qualitatively the cutting forces and 
vibrations of the cutting tool depend on the angle of inclination of the surface. 

An analysis of the experimental data made it possible to identify a trend towards a decrease in the amplitude 
parameters of roughness by 10—20 % with an increase in the angle of inclination to 30-40°. The quality of the treated 
surface was reduced due to the elastic retraction of the spherocylindrical tool and an increase in the amplitude of 
vibrations of the technological system. 

It was established that the forces grew in direction a,(X) and fz(Y) at y > 40°. This can be explained by the 
distribution of the cutting force components along the cutting edge, which depends on the inclination of the surface. 

The solutions to the problems of monitoring and analysis of roughness parameters obtained by the authors due to 
acoustic diagnostics can significantly reduce the volume of experimental studies. This approach makes it possible to 
significantly reduce the cost of material resources and time for mechanical processing of products. 

The experimental results are of practical interest and make the information about the formation of microrelief, 
including the roughness parameters under mechanical cutting, more accurate. A more detailed study of the proposed 


method of acoustic diagnostics is needed. 
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Abstract 
Introduction. Some of the major factors in the occurrence of vibration of units, causing fatigue failure of the housings 


of elements, pipelines, and failure of pump elements, are pulsations of the working medium in the hydraulic systems of 
machine tools, fuel feed systems of aircraft engines and liquid-propellant engine supply. This study aimed at the 
implementation of a method for determining the dynamic characteristics of a volumetric pump using special bench 
systems, and the comparison of the calculation results to the experimental data. The stages of calculating the dynamic 
characteristics of a volumetric pump were described, taking into account the pre-developed special bench systems on 
the example of an external gear pump with a capacity of 14 cm*/rev. The implementation of V. P. Shorin's load 
variation method using special bench systems developed by the authors with predetermined dynamic characteristics was 
shown. The main stages of the methodology for determining the dynamic characteristics of a gear pump were described. 
Materials and Methods. Methods of spectral analysis of pulsating pressure were applied in the work. Pulsations 
of fluid flow at the pump outlet were determined using the impedance method, the method of load variation, and 
special bench systems. 

Results. The paper implemented a technique for determining the dynamic characteristics of a gear pump in the drive 
shaft speed range of 500-2500 rpm for four harmonic components of the vibration spectrum in a wide range of dynamic 
loads (from inertial to capacitive nature). The bench systems yielding the calculation of the dynamic characteristics of 
the pump with a minimum error based on the condition of matching the dynamic load and the source of vibrations were 
analyzed. The developed approach to the evaluation of the dynamic characteristics of the pump was verified through 
comparing the calculated and experimental data of pressure pulsations in the bench systems with choke, cavity and an 
extended pipeline at the pump outlet. 

Discussion and Conclusions. The method for determining the dynamic characteristics of a volumetric pump was 
implemented using special bench systems developed by the authors. The research results show that the gear pump under 
study can be considered as an independent source of flow fluctuations, for which the deviation of its own dynamic 


characteristics from the average values does not exceed 10% for the first harmonic component. 

Keywords: volumetric pump, load variation method, bench systems, source of flow fluctuations, pump dynamic 
characteristics, independent source of fluctuations, flow pulsations, impedance. 

Acknowledgements: the authors would like to thank the reviewers for the comments, which contributed to quality of 
the article — the research of the second author was carried out on theme no. FSSS-2020-0015 within the frame of the 
government task of the Ministry of Education and Science of the Russian Federation. 

For citation: V. I. Sanchugov, P. D. Rekadze. Determination of the dynamic characteristics of a gear pump by the load 
variation method using special bench systems. Advanced Engineering Research, 2022, vol. 22, no. 2, pp. 130-141. 
https://doi.org/10.23947/2687-1653-2022-22-2-130-141 


© Sanchugov V. I, Rekadze P. D., 2022 


Sanchugov V. I., et al. Determination of the dynamic characteristics of a gear pump by the load variation method using special bench systems 


Introduction. The main disadvantages of pumps are the acoustic noise that accompanies their operation, the 
vibration of the body, connected pipelines, and the resulting structural damage. 

The traditional way of improving pump designs, used in mechanical engineering and consisting in smoothing out the 
uneven flow of the liquid, is hampered by the lack of means of direct measurement of fluid flow at the pump outlet. 
This limits any experimental evaluation of the design efficiency and finishing work on pumps. In addition, this limits 
the analytical description of the pump as a source of vibrations required by hydraulic system developers to provide the 
target dynamic quality of the systems by themselves. 

For this reason, the research on the modeling of dynamic processes in hydromechanical systems [1,2] and the 
application of new approaches to solving problems of hydraulic system dynamics [3, 4] based on physical principles are 
pressing issues. 

Materials and Methods. The implementation of V. P. Shorin's method of load variation [4] is presented. It consists 
in determining pressure pulsations behind a pump interacting with several bench systems, whose dynamic 
characteristics are calculated in advance. Further, the processed frequency-response functions of the pressure using the 
impedance of the bench systems are recalculated into the pump's own dynamic characteristics — pulsation performance 
and impedance. The development of the approach provides that special bench systems are proposed. They take into 
account the connecting fittings, adapters and internal channels of the pump, which makes it possible to more accurately 
determine the dynamic characteristics of pumps, as well as to verify the applicability of the pump model used in the 
form of an equivalent source of vibrations. 

Research Results 

Method for determining the dynamic characteristics of a gear pump 

Samara University has developed a computational and experimental method for determining the dynamic 
characteristics of a gear pump (pulsation performance and impedance), which is based on the research of L. Thévenin, 
E. Norton and V. P. Shorin [5, 6]. The developed methodology includes the following main stages: creation of models 
and calculation of dynamic characteristics of special bench systems; experimental determination of pressure pulsations 
at the pump outlet; calculation of the spectrum of excited vibrations in various bench systems; determination of 
approximating dependences for individual harmonic components; calculation of the pump own dynamic characteristics 
(pulsation performance and impedance); evaluation of the independence (stability) of the dynamic characteristics of the 
pump from the characteristics of bench systems; verification of the applicability of the pump model used as an 
equivalent source of flow (or pressure) fluctuations. The following is an example of the implementation of the proposed 
technique using an external gear pump with a specific capacity of 14 cm*/rev. 

Models and calculation of dynamic characteristics of special bench systems 

Basic options of bench systems at the liquid outflow during pump tests were considered. In concentrated parameters, 
the bench system can implement active flow-dependent, inertial, and elastic loads. In distributed parameters, the bench 
system can be provided in the form of an extended cylindrical pipeline implementing an active frequency-independent 
non-reflective load. 

Pump models in the form of equivalent source of flow fluctuations (ESFF) or pressure (ESPF) are determined using 


constant flow source Q,, an ideal source of pressure fluctuations p (or flow q ) and internal impedance Z, [4]. Figure 1 


shows pump models in the form of an equivalent source of pressure fluctuations (and flow) with an attached system. 


bench 
pump | system pump | system 


a) b) 


Fig. 1. Pump models in the form of an equivalent source of: a — pressure fluctuations; b — flow with attached system 


Q,,P, — average flow and pressure; p,,,g,, — pump's own pulsating performance in terms of pressure and flow, 


Z.,,»Zem.c. —~ impedance of pump and bench system (adapted from [4]) 


ur 
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In this connection, the task of selecting a bench system is to provide such system parameters at which the best 
conditions for high measurement accuracy are realized. The task of selecting the parameters of the attached system 
should be carried out when the impedance changes in the range from -00 to +00, 

Figure 2 shows schematic diagrams of the bench systems. 
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Fig. 2. Schematic diagrams of the bench systems for determining own dynamic characteristics of pumps: a — with active throttle 
load; b — with active non-reflective load; c — with capacitive load: 1 — test pump; 2 — load choke; 3 — flow meter; 
4 — pressure sensor; 5, 6 — connecting fittings; 7 — main line; 8 — extended pipeline; 9 — cavity 


Computational models of impedances of special bench systems: 
— with active throttle load (throttle): 


ons, LE (1) 
. 1 Jet anes, 
Zo (Q,, f) = |(Z,, (Q, ) +ReZ,,, (Q, yy [ann dl e  Zin(Qo)+ReZnp( Qo) 


mp 


— with active distributed load (extended pipeline): 


pa (2) 


— with capacitive load (cavity): 


‘ (3) 
2 | (2af,) LC-1 
Zone (Qos fi) = q|( REZ ses (Qo) + REZyys (Qn) + REZ (Op) + Zip (Qa)) a ; 
k 
wr (27h; y Lc-1 
27 ficC(REZ pay, +RCZ nypa (Sic + REZ nc. (Fi }+Zop (Qo )) 


The formulas indicate: Z,, (Q,) — impedance of the throttle; ReZ,,(Qo) — real part of the impedance of the 
pipeline, including the connecting line, fittings, internal channels in the pump and _ throttle; 
ReZ nam (Qo), REZryp6(Qo), ReZu.c.(Qo) — real part of the impedance of the fittings and the internal channel in the 
pump and cavity with laminar, turbulent fluid flow behind the pump and local resistances, respectively; f, — 
frequency k — th of the harmonic component of the vibration spectrum behind the pump;  — density of the working 
environment; /, — pipeline length, including the length of main line 6, connecting fittings 5 and 7 and internal 


channels of the pump and throttle (not shown in Fig. 7); S,,, — cross-sectional area of the pipeline; L — “inertia” of 
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the connecting fittings and the internal channel in the pump and cavity; C — elasticity of the cavity; k — number of 


the harmonic component of the vibration spectrum; j — imaginary unit (j = |. 


Using a bench system with a throttle at the pump outlet, three inertially active dynamic loads were implemented, in 
which the active resistance was changed through changing the overlap area of throttle Fo: 

— for bench system no. 1, overlap area Fo, =9.290x 10° m’, 

— for bench system no. 2, overlap area Fy:=7.299x10° m’, 

— for bench system no. 3, overlap area Fo3=2.247x10° m’. 


Using a bench system with a cavity at the outlet of the pump with a volume of 2 liters, an inertial elastic dynamic 


load was realized. 


An active frequency-independent load was implemented using a bench system with an extended pipeline 


length of 106 m. 


The dependences of the impedance (module and phase) on the frequency were calculated according to the models of 
bench systems using formulas (1)—(3). Table 1 for the Ist harmonic shows the regression functions of the impedance 
module of the bench system Z;(f) and phases B;(f) for the analyzed bench systems. 


Table | 


Regression functions of dynamic characteristics of bench systems according to the 1st harmonic 


No. of stand Bench system impedance 
De Hensnaytomee module Bench system impedance phase §; (f), 
according to Z,(f) ke rad 
Figure 2 mas? 
a—no.1 ainciioat 1.412-10’f + 2.201-10" 7.066: 10™f + 0.199 
a—no.2 a 1.610°10'f + 4.25810" 6.906- 10°F +0.09039 
a—no.3 4.980-10'f + 1.503-10” 2.723-10°f + 0.251 
.308-10"f? — 1.416-10’ 4.234-10° f +1.495-10° 
caw With cavity at pump 3.308:10°f — 1 ale 10°'f+ Grete - f ; 
outlet 1.873-10 1.430-10° f +2.930-10 
With extended 
c —no.5 pipeline at pump 1.897-10'° 0 
outlet 


Experimental determination of pressure pulsation at the pump outlet 


The background information for determining the dynamic characteristics of the pump was the measurement of 
pressure at the pump outlet. An example of an oscillogram and a pressure pulsation spectrum at the pump outlet in the 
bench system with an active distributed load is shown in Figure 3. Operating mode of the pump in this case is as 
follows: average pressure PR, = 23.8 MPa, drive shaft speed n=1000rpm. 


24.30 ppssseeeseee ae Pmax=24.26 Mpa|------------== prevessssseeeed persssseeeeay} 
24.13 | : | . | : 
23.78 | 


Pressure, MPa 


OA AS isch feasted: decal Pesaat df cactac i a ¢ faces | weed + Soe Ay Lp 
Q395 Miscel| snceod ue seeehelhedteselfeed a io (eee tp 


P3084 heck eee ee one ee ee A Meee 


Pmin=23.05 Mpa ‘ 


T s T Y T Y T —sa 
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Rotor rotation time, s 
t T T T T T 1 
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Rotor rotation angle, deg 


a) 


Machine building and machine science 


133 


http://vestnik-donstu.ru 


Advanced Engineering Research 2022. V. 22, no. 2. P. 130-141. ISSN 2687-1653 


Amplitude of pressure pulsations, MPa 


167 334 501 668 835 1002 1169 1336 1503 1670 1837 
Frequency of oscillation spectrum components, Hz 


b) 
Fig. 3. Pressure at the pump outlet: a — oscillogram; b — spectrogram 


The amplitude-frequency response of pressure pulsations at the pump outlet in the speed range of the pump drive 
shaft n=500—1250 rpm is shown in Figure 4. The regression functions built in the Microsoft Excel environment using 
the Trendline tool are plotted with a solid line. 


® | harmonic 
® 2 harmonic 
@ 3 harmonic 


e@ 4 harmonic 


Amplitude of pressure pulsations, MPa 


0 83 166 249 332 415 498 581 664 747 830 


Frequency of vibration spectrum components, Hz 
———— 3 


0 500 1000 1500 


Frequency of power shaft rotation, rpm 


i: -_ 
0 73 146 21,6 


Average pump fluid flow, I/min 


Fig. 4. Amplitude-frequency response of pressure pulsations at the pump outlet by individual harmonic components 
(dots indicate data on: ® — first, ® — second, ® — third, @ — fourth harmonics, lines indicate regression functions) 


The analysis of the oscillogram, spectrogram and amplitude-frequency response of pressure pulsations showed the 
following: the type of process was polyharmonic steady; the number of recorded harmonics was 8, analyzed! — 4; the 
amplitude range of pressure fluctuations was 0.28—0.52 MPa; the drive shaft frequency range was 500-1250 rpm, the 
analyzed harmonic components of the spectrum 83-833 Hz. 

Fourier series expansion was used for harmonic analysis of periodic signals [7]. To decompose into a Fourier series, 
sections of the pressure oscillogram measured relative to the reference signal at the rotor frequency were used, i.e., one 
section of the oscillogram was extracted for a full rotation of the rotor. The reference signal was recorded at a fixed gear 
position, set at the beginning of recording the pressure waveform for five bench systems. In this way, the amplitude- 
frequency (FRF) A;,(w) and phase-frequency (PFC)@;,(w) characteristics were determined. 


' Tt is assumed in the work that the extreme highest analyzed harmonic in amplitude should not be more than 5 % of the amplitude of the first 
harmonic. 
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Calculation of the spectrum of excited vibrations in various bench systems and determination of approximating for 
individual harmonic components of pressure fluctuations 

Based on the experimental determination of pressure in five bench systems: with a throttle (no.1—3), with a cavity 
(no.4) and an extended pipeline (no.5) at the pump outlet, the spectra of excited vibrations were calculated. Table 2 
shows the regression functions of the amplitudes A,(f) and the initial phases of pressure pulsations for the 1st harmonic. 


Table 2 
Regression functions of amplitudes and initial phases of pressure fluctuations 
No. of stand sas 
system Pressure pulsation amplitude Initial phase ph Frequency 
. Bench system type kg pressure pulsations 
according to Api &), | tad range 
Figure 2 be 
—no.1 1.003-10°f 0.543 
a With throttle at : 
a—no.2 1.243-10°f —0.280 
pump outlet x 
a—no.3 6,129-10°f —5.941 
With cavity at pump 1.343-10°F — 4.298-107f + 83-417 
poe outlet 3.913-107 a 
With extended 
c —no.5 pipeline at pump 2.694: 10°F —0.658 
outlet 


Calculation of the dynamic characteristics of the pump 
To calculate the dynamic characteristics of the pump (impedance and pulsation performance of the source), it is 
required to conduct an experiment with at least two dynamic loads”. Therein, according to the model of an equivalent 


source of flow fluctuations refined by the authors using analytical dependences p, ( f ) and Z, ( f ) from Tables 1 and 2, 


the pump impedance is calculated from the formula: 


jon 
P, 
L =) 4 (4) 
PL 
P, Z, 
where Z, — pump impedance; ( Z, = [Zee ) — pump impedance module; (|Z,|=./ReZ,,° + ImZ,," ), ReZ, ,lmZ,, — 
real and imaginary part of pump impedance; @, — pump impedance argument; Z,,Z, — input impedance of dynamic 


loads no. 1-2; (Z, =|Z,|e", Z, =|Z, |e"), 2,8, — impedance arguments of dynamic loads no. 1-2; p,,p, — 


pressure pulsations behind the pump in bench system no. 1-2; ( p, = A,,e"", p, = A,,e’” ), A,,,A,. — amplitudes of 


flow pulsations behind the pump in bench system no. 1-2; g,,g, — initial phases of pressure pulsations of dynamic 
loads no. 1-2. 

Then, using the analytical dependences from Tables 1 and 2, the pulsation capacity of the pump was determined 
according to the model of an equivalent source of vibrations refined by the authors (variable component of volumetric 
flow q, and pressure p, ): 


ma 14 
P2 Z, Z, 
uo ? Py = (5) 
Z, l P2 PZ, 
P P, 2, 


? For a sufficient sample size, it was required to use at least four dynamic loads Z; , that form six non-repeating (including i # j) combinations: Z,, , 


Zi3» 2g Za3» Loy, Lrg 
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where q,,p,, — variable component of the volume flow (gq, = Ae") and pressure ( p, =A,,e""") of the pump; 


A,,»A,,— pump flow and pressure pulsation amplitude; g,,~, — argument of the variable component of the volume 
flow and pressure of the pump. 
Of the calculated characteristics according to ESFF (q, (f)) and ESPF ( p,(f)) models, using formula (5), only 


one will be stable, since the equivalent oscillation source has an independent dynamic characteristic only for one of the 
parameters (pulsating flow or pressure). 

The calculation of impedance Z, and the variable component of the volume flow g, of the pump was carried out 
according to formulas (4)—(5) for five dynamic loads, i.e., using ten different combinations of these loads: “1_2”, “1 3”, 
“237,14, 24, 34 “15”, 25”, “35”, “45” (Table 1). 


For further calculation and analysis, a characteristic of the pump flow pulsation amplitudes A, was obtained from 


ten flow pulsation amplitude curves. The values of the amplitudes A, (at each frequency) were calculated using the 
formula from source [8]: 

Aw = Au £ BA as 

where Au — average amplitude of pulsating flow; AA,, — confidence interval of the amplitude of pump flow 


pulsations (AA, = aes, o — standard deviation of the amplitude of the pulsating flow rate from the average value; 


qu in 


m — number of dynamic load combinations (m=10); t— Student's coefficient (t = 2.262, was calculated for the 
confidence level 0.95). 


The confidence interval of amplitudes AA, was formed from errors (instrument and method) and the degree of 


dependence of the dynamic characteristics of the pump A, ( f ) on the bench system. The instrument error consisted of 


the error of determining the pressure and phases, the method error included the error of the unaccounted influence of the 
dynamic properties of the drive mechanical system. 

The phases of flow pulsations gy, and their deviations were calculated in a similar way. When calculating the 
phases, combinations “1 2” and “1 5” were excluded due to the degeneration of the variable component of the volume 
flow of the pump g, when using close values of characteristics at the same frequency: 9, ~ @, and f, = f, . 


The results of calculating the intrinsic dynamic characteristic of the pump in the form of pulsation performance and 


impedance are shown in Figure 5. Phases, for convenience, are presented in degrees. 
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d) 
Fig. 5. Dynamic characteristics of the pump: a — flow pulsation amplitude; b — flow pulsation phase; c — impedance modulus; 
d — impedance phase (solid lines indicate data on: eum — first, @======—" — second, 
eum — third, emmmmme — forth harmonics) 
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The amplitude of the flow pulsations at the 1st harmonic has a monotonically increasing character. With an increase 
in the harmonic number, the slope of the flow pulsation amplitude curve decreases to an almost constant value to the 4th 


harmonic. The values of amplitudes A, at the higher harmonics are by an order of magnitude less than at the Ist 


harmonic. 
The pump impedance module at the Ist harmonic has a monotonically increasing character. 137 
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To evaluate the realized dynamic loads, the criterion of matching the impedance module \Z,,| of the oscillation 


was used [99, 10]. The impedance modules should 


cm.c 


source and the input impedance module of the dynamic load |Z 


be comparable (|Z, ), so that the quantitative values of the dynamic characteristics of the pump were 


=|Z 


cm.c 


u 


u 


was used. At 


determined with minimal error. For comparison, ratio ~1, the dynamic load was considered 


cm.c. cm.c. 


consistent with the source. 
Due to the difference in the curve shapes of the pump impedance modules and dynamic loads in the frequency range 


u 


of 83-417 Hz, it was rational to use average ratio , that provided considering the contribution of most of 


cm.c. 
min **Fmax 


u 


the frequency range, i.e., ratio according to the Ist harmonic would be: 


cmc. 


Finin -+-Simax 

— for active throttle load no.1 — 1.16; 

— for active throttle load no.2 — 0.8; 

— for active throttle load no.3 — 2.65; 

— for capacitive load no.4 — 4.64; 

— for active distributed load no.5 — 0.35. 

That is, the dynamic loads of bench systems no. 1, 2 were the most consistent with the source of vibrations. Thus, 
the dynamic characteristic of the pump, calculated using dynamic loads of bench systems no.1, 2 (“1_3”, “2.3”, “1_4”, 
“2 4”, “1_ 5”, “2 5”) was determined with a minimum error. 

Evaluation of the stability of the dynamic characteristics of the pump from the characteristics of bench systems 


Stability of the dynamic characteristics A, and 9, (IZ, and @, ) of the pump was evaluated using mathematical 


Statistics tools [8]. Thus, the stability of the characteristic A,, was estimated through calculating the ratio of the 


confidence interval |AA at this frequency* The 


qu 


at frequency f to the average value of the characteristic [Anat 


resulting relative deviation A, should not exceed the specified deviation A,,,, : 


ae MA, on 
q | | — dong? 
quMe 
where [A uae — average median value of the amplitude of pump flow pulsations; A,,,, — specified relative deviation 


of the pump flow rate pulsation amplitude from the average value. 


With relative deviation A, 0 (when |A4,, —> 0), the stability of characteristic A,, was maximum. Based on the 


experience of other researchers, it is preferable to select value A,,,, within 5-30 % [11-17]. 


For a single-digit evaluation of the dynamic characteristics of the pump, it is rational to use the value of the relative 
deviation. The calculation of such a value is possible through estimating the average median value of relative deviations 


AA 
in the frequency range of each harmonic ( fyi,---fmue): Mgve =—— . Value A,,,, is efficiently used, if 
PN fin Sma 
MA, . . : . . 
< 0.3, i.e., when the confidence interval AA,,, does not exceed 30% of the average amplitude of pulsating 
A 
Fin Sina 


qu * 
The results of assessing the stability of the dynamic characteristics of the pump from the characteristics of the bench 
systems’ were: according to the Ist harmonic — 2 %, according to the 2nd — 8 %, according to the 3rd — 15 %, 
according to the 4th — 48 %. The results of the calculation of the dynamic characteristics of the pump showed that for 


> Stability of other intrinsic characteristics of the pump is evaluated similarly. 
“Tt is required to consider the method error of the approach and the instrument error of determining the pressure. 
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four harmonics in the frequency range of 500—2500 rpm, the gear pump under study should be considered according to 
the model of an equivalent source of flow fluctuations. 

Checking the applicability of the pump model used as an equivalent source of flow fluctuations 

To test the developed approach for evaluating the dynamic characteristics of the pump, a calculated determination of 
pressure pulsations in bench systems with a throttle, a cavity and an extended pipeline at the pump outlet was 
performed. Using the found dynamic characteristics of the pump, the amplitudes of pressure pulsations in the analyzed 
bench systems were calculated, which were compared to the experimental ones. The results of the calculation of 


amplitudes A, ( f ) in a bench system with an extended pipeline are presented below. The amplitudes of pressure 


pulsations p, were calculated in accordance with the ESFF model from the formula: 


qi . Z, . Z; (6) 
| aaa 
Z,+Z, 
The results of calculating the amplitudes of pressure pulsations behind the pump are shown in Figure 6. The dotted 
line indicates the amplitudes obtained experimentally, and the solid line — calculated from formula (6). 
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Fig. 6. Pressure frequency response for the bench system with extended pipeline length /=106 m (solid lines indicate calculated data 


On: Gomme — first, Ge — second, eum — third, em — forth harmonics; 


dotted lines indicate experimental data) 


Convergence of the results of the verification calculation was evaluated according to the formula of deviation of the 
difference between the experimental and calculated values from the experimental data. The results of the test calculation 
for the bench systems with a throttle, a cavity and an extended pipeline at the pump outlet showed that according to the 
1st harmonic, the relative deviation from the experimental data for the considered bench systems did not exceed 10 %, 
the minimum relative deviation was characteristic of a bench system with an extended pipeline at the pump outlet — 
2%. The amplitudes of the pressure pulsations of the 2nd—4th harmonics were within the corridor calculated from the 


formula A, =A,+A, (where A, — the average amplitude of the pulsating pressure). 


Discussion and Conclusions. Methodology for calculating dynamic characteristics of the pump in the form of an 
equivalent source of fluid flow fluctuations was implemented on the example of a gear pump with a capacity of 14 
cm’/rev using special bench systems developed by the authors and an updated model of impedance and pulsation flow 
rate. It was shown that the investigated gear pump of external gearing was rationally considered according to the model 
of an equivalent source of flow fluctuations. The proposed approach was verified in the form of calculation of pressure 
pulsation amplitudes in the bench systems and comparison to experimental values. The deviation of the calculated 
amplitudes of pressure pulsations from the experimental values did not exceed 10%. It is worth noting that the impact 
on the dynamic characteristics of the gear pump of the drive mechanical system and the wear of the gear pump parts 
remained outside the scope of the review. Based on the work carried out, it was possible to formulate directions for the 
development of the proposed approach: 

— assessment of the dynamic characteristics of the gear pump during the entire service life; 
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— development of a four-pole model of a gear pump as an equivalent source of vibrations, taking into account the 
mechanical drive and hydraulic output systems; 

— development of a gear pump model as an equivalent source of vibrations, taking into account its design features, 
the displacement process, reverse water hammer in the chamber between teeth, the work of the locked volume, and 
leakage from the injection area to the suction area. 

Here, it is required to expand the nomenclature of the tested volumetric pumps for the entire class of volumetric 
displacement machines to confirm the operability of the proposed approach. 
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Abstract 

Introduction. A gravity conduit with a control valve and a pressure regulator in the gravity water-supply pipe is 
considered under conditions when the hydraulic regime of the water-supply system is not disturbed. In relation to such a 
system, the problems of converting the energy of local artificial hydraulic resistance into electricity are investigated. 
Materials and Methods. Literature that highlights the possibility of using microturbines for power generation in water 
systems was studied. The actual values noted by the continuous pressure recorder (logger) for 12 hours were presented, 
as well as the water consumption by the turbine at a given unit section (the average for a year), and the pressure 
differential. It was noted that the use of small hydroelectric power plants in water supply systems significantly reduced 
the cost of their operation. The indicators of water consumption in the hydroturbine unit during the year, broken down 
by months, were given. The maximum power at the turbine inlet was calculated. The principles of selecting the type of 
hydraulic turbine were described. The average efficiency values for different elements of the hydroturbine plant, the 
average parameters of the power of the small hydro and the corresponding indicators of the average monthly electricity 
generation were indicated. 

Results. Equipping the units under study with specially designed turbines can enable to obtain electrical energy through 
converting artificially created by local resistance and extinguished mechanical energy. It is possible to apply the 
approaches described in this article when replacing many of the pressure control units of the Yerevan City Network 
system. The productivity of a small hydro power plant was predicted, as well as the terms of its construction and 
operation — for 2 years and 30 years, respectively. The construction and maintenance costs were calculated in advance. 
The expected data on income, expenses, and net profit are given. It is indicated that it will potentially be close to 6 
million drams per year. The analysis of the data allowed us to conclude that the internal rate of return will be at the level 
of 10.4%, and the payback period is 9 years. 

Discussion and Conclusions. We recommend replacing the regulator with a hydroturbine having the same hydraulic 
resistance and automatic flow control in the system. The conclusions were confirmed by the energy and economic 
indicators of the hydroturbine plant located on the section of the Arzakan — Yerevan main water pipeline. 
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Introduction. Traditional non-renewable primary energy sources, which are very limited in resources, continue to 
play a dominant role in the Armenian economy. The active use of fossil energy sources has caused environmental 
problems such as air and soil pollution, water scarcity, and ecosystem degradation. Population growth, continuous 
household and industrial energy consumption should be taken into account. All this suggests that further quantitative 
(extensive) development of the energy industry is inexpedient. The use of renewable energy sources to replenish the 
overall energy balance of the country is an imperative of the time. This approach makes it possible to provide 
sustainable development while preserving the environment. Water and energy resources are essential for human 
life and are subject to regular economic, technological, demographic, and social impacts [1]. According to some 
estimates, 2—3 % of the world's electricity consumption comes from water systems (WS) [2], and pumping units are 
used in 80—90 % of cases. The costs of their maintenance are among the main operating costs of WS [2]. The concept of 
rational and efficient use of water resources and electricity is of strategic importance for sustainable development and 
climate change mitigation. Rational use of water and electricity is hampered by weak infrastructure and outdated 
operating procedures. This is especially true for developing countries. Moreover, in line with the Millennium 
Development Goals', there is a need for alternative solutions: the goal is to halve the number of people without 
sustainable access to clean water and basic sanitation. As renewable energy sources, it is proposed to put into operation 
small hydroelectric power plants (SHPP) in sections of the water supply network with local hydraulic resistance, subject 
to semi-open valves or pressure regulators. This is cost-efficient for the production of electricity and water supply, and 
illustrates perfectly the rational use of water resources. In the past decades, hydroelectric power was one of the main 
sources of energy from water. Today, the water sector is considered a direct consumer of electricity, and this affects the 
distribution of water resources’. Active energy consumption directly affects the state of water resources in the world, 
which in turn is the cause of climate change. Rational use of energy should contribute to sustainable development 
through the competent distribution of energy resources at all stages of transformation. For a comprehensive assessment 
of the problem, an example of the construction of the SHPP on the Arzakan —Yerevan water pipeline is given. 

Materials and Methods. The Arzakan-Yerevan water pipeline supplies water to the administrative districts of 
Yerevan (Arabkir, Malatia-Sebastia, Achapnyak, Erebuni) and adjacent settlements (Zovuni, Kanakeravan, Nor-Achin, 
Nor-Gekhi, Yeghvard, etc.). The water pipeline originates from the Arzakan spring, and on the territory of the Getamej 
administrative community, it is divided into two water pipelines with different pressure conditions. The construction of 
the SHPP on one of the branches of this water pipeline feeding the north-western and western districts of Yerevan is 


being considered. The company operating on the dividing branch of the water pipeline installed a semi-open valve to 


' Millennium Development Goals. United Nations URL: https://www.un.org/development/desa/ru/millennium-development-goals.html (accessed: 
05.06.2022) (in Russ.) 

’Mao Yushi, Sheng Hong, Yang Fuqiang. The true cost of coal. Understanding China's Energy Landscape. 2008. URL: 
https://www.understandchinaenergy.org/greenpeace-chinathe-china-sustainable-energy-programwwft-china-main-authors-mao-yushi-sheng-hong- 
yang-fugiang-2010-the-true-cost-of-coal-full-report/ (accessed: 05.06.2022). 
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adjust the water flow in the system. This provides the pressure in the water supply for settlements located above the 


main highway (adjacent to the city of Nor-Achin). When the valve is semi-open, a local hydraulic resistance of about 


13 m of water column (w.c.) is created (before that, the pressure in the water pipe was about 33 m w.c.). 


We propose to install a SHPP instead of a semi-open valve, without disturbing the normal hydraulic mode of 


operation of the water supply system. Figure 1 shows a plan of the area in the Google Earth program, and Table 1 


summarizes the actual values recorded by the continuous pressure recording device (logger) for 12 hours. 


Data from pressure recorders installed in the study area 


Fig. 1. Map of the area in Google Earth 


Table 1 


Research Date ime Registered pressure (m H20) 
Before semi-open valve After semi-open valve 
15.09.2020 20:00 32.8 13.0 
15.09.2020 21:00 33.2 12.6 
15.09.2020 22:00 33.3 13.0 
15.09.2020 23:00 33.3 13.0 
16.09.2020 00:00 33.5 12.8 
16.09.2020 01:00 33.8 12.4 
16.09.2020 02:00 33.9 13.2 
16.09.2020 03:00 33.9 13.1 
16.09.2020 04:00 33.9 13.6 
16.09.2020 05:00 33.8 13.3 
16.09.2020 06:00 33.8 13.2 
16.09.2020 07:00 33.7 12.6 
16.09.2020 08:00 33.4 13.0 


The proposal aims at providing the transformation of hydraulic energy artificially created in the local resistance unit 


into electrical energy. Favorable conditions are created at the construction site of a small hydroelectric power plant (on 


a branch of the Arzakan-Yerevan main water pipeline) for the integrated use of water and energy potential. The flow 
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rate of the turbine in this section of the assembly is 1.15 m/s (average during the year). The control device on the pipe 
provides a given flow rate. As a result, pressure drops from 20.3 m w.c. to 21.4 m w.c. (Table 1). The expediency of 
the practical use of the so-called “hidden” energy generated in the sections of main water pipelines in the Republic of 
Armenia should be noted. The hydropower potential of water supply systems has long been known, but not properly 
explored. The literature describes in detail the possibilities of using microturbines to generate electricity in water supply 
systems [3-5]. Gravity systems in areas with high slopes form significant pressures in water supply and distribution 
networks, creating the prerequisites for generating hydroelectric power. Turbines installed in water distribution 
networks can also be used as pressure control systems. To do this, pressure relief valves (PRV), which control water 
losses and leaks, are used [6]. They dissipate energy, thus reducing the pressure in the system. As for water turbines, 
they can convert excess pressure into electricity®. Major advantages of hydraulic energy recovery in water pipes, 
according to F. Vieira and J. S. Ramos [3, 7], are improving the energy efficiency of the system through the use of local 
sources, and reducing dependence on external (grid) energy. In addition, hydraulic energy recovery contributes to an 
overall reduction in operating costs. F. Vieira and J. S. Ramos [3, 7] emphasize in their research that the use of small 
hydroelectric power plants in water supply systems reduces significantly the cost of their operation. This is due to the 
fact that the proposed solution enables to replace the control valve with a hydraulic turbine equal in hydraulic 
resistance. Table 2 shows the water consumption in the system (monthly averages are indicated). 

Table 2 


Water consumption in the hydroturbine unit 


Month Consumption, m°/s 
January 1.20 July 1.30 
February 1.10 August 1.30 
March 1.10 September 1.05 
April 1.00 October 1.05 
May 1.10 November 1.10 
June 1.30 December 1.20 


Note that the maximum head is Hmax = 21.4 m, the minimum — Hmin = 20.3 m. 
Average annual consumption, according to Table 2: Q = 1.15 m’/s. 
With a design head equal to H = 20 m, the flow power at the turbine inlet [12]: 


N= pqQH = 1000 - 9.81 - 1.15 - 20 = 225632 Br 225 kW. 


Principles for selecting a hydro turbine. With a low head and a higher flow rate (for SHPP), a turbine with a high 
specific speed should be selected. Of the reaction turbines, an axial jet turbine will do. Of the low-power active turbines, 
Banki turbine should be chosen to pass a relatively large flow rate. However, in this case, it is unacceptable to use an 
active turbine, since when it is filled, excess pressure is formed at the outlet, and Banki turbine operates as a very poor 


jet turbine with very low energy performance. 


> The Millennium Development Goals Report 2011. United Nations. New York, 2011. URL: 
https://www.un.org/millenniumgoals/pdf/(2011_E)%20MDG%20Report%202011 Book%20LR.pdf (accessed: 05.06.2022). 
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It is recommended to install the hydro turbine unit at the valve on the bypass line. Inlet and outlet valves will 
provide uninterrupted operation of the line even in the event of a station shutdown. The operation of the SHPP will be 
monitored by pressure sensors installed at the inlet and outlet of the turbine. This maintains the required pressure on the 
line*. When installing equipment on drinking water mains, special requirements should be taken into account [8-10]. 
Specifically, contamination of the water with lubricating oils and other materials” ° must be avoided. In general, the 
distributed water should meet the quantitative and qualitative standards of drinking, domestic, and industrial 
consumption [11]. 

Energy indicators of SHPP. The following average efficiency values are applied for different elements of the 
hydro turbine installation: n,, =0.86 , ,,, =0.94. 

Turbine shaft power: N, = n, N = 0.86 225 = 193.5 kW, and at the generator output: Nrex = Nren Nz = 0.94 » 193.5 = 
184 kW. 


Table 3 presents the averaged parameters of the power of the SHPP and the corresponding indicators of the average 


monthly electricity generation. 


Table 3 
Average monthly capacity and output of SHPP 
Month Power, kW Output, kWh 

January 190 141,360 
February 174 116,930 
March 174 129,456 
April 158 113,760 
May 174 113,456 
June 206 148,320 
July 206 153,264 
August 206 153,264 
September 166 119,520 
October 166 123,504 
November 174 125,280 
December 190 141,360 
Total 1,573,474 

It should be noted that the use of a SHPP jet turbine at a large negative suction height (A, =—13m), in our 


opinion, is proposed for the first time in world practice. In this case, on the one hand, the occurrence of cavitation is 


“Markaryan AYa, Tokmadzhyan VO. Regulirovanie proizvoditel'nosti nasosov s tsel'yu predotvrashcheniya kavitatsionnykh yavlenii. In: Proc. 7th 
Int. Congress “Water: Ecology and Technology”. Moscow: EHK VATEHK-2006. Part 1. P. 566. (In Russ.) 

°Sanitarakan Kanonner yev Normer N2-III-A2-1, Khmelu jur: Jramatakararman kentronats’vats hamakargeri jri vorakin nerkayats’vogh higiyenik 
pahanjner, Voraki hskoghut’yun, Yerevan, 2002t, ej 11. (In Armenian) 

° Barry JA. Watergy: Energy and Water Efficiency in Municipal Water Supply and Wastewater Treatment Cost-Effective Savings of Water and 
Energy. The Alliance to Save Energy. Washington, 2007. 44 p. 
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excluded [9, 12], which, undoubtedly, is regarded positively. On the other hand, it is not known how the high 
pressure at the outlet of the turbine affects the efficiency, therefore, in our calculations, a slightly underestimated 
indicator is used. 

Research Results. The electricity produced by the SHPP (Table 3) is planned to supply to the general energy 
system of the Republic of Armenia on the terms of guaranteed purchase at prices approved by the Commission on 
Regulation of Public Services of Armenia (currently, the tariff is 10.579 drams per kWh). 

According to Table 3, during the year, the SHPP will produce at least 1.5 million kWh. According to preliminary 
calculations, the construction period of the SHPP is 24 months, and the service life is 30 years. The cost of construction 


and design of the SHPP will amount to 60.4 million drams (Table 4). 


Table 4 
SHPP construction costs 
Works Cost, AMD mln 
Pipeline construction / maintenance 2.2 
Construction of the HPP building 14.5 
Cost of installing hydraulic units (hydro turbines, generators) 18.0 
Transformer substation 10.7 
Construction of high voltage power lines 3.0 
Design and calculation work 7.0 
Other costs 5.0 
Total 60.4 
Anticipated revenues, operating costs, and profit are shown in Table 5. 
Table 5 
Forecast of annual financial indicators 
Indicator Amount, AMD min 
Gross income 15.9 
Salary, operating expenses, current repairs, and elimination of accidents 7.0 
Duties and other obligatory payments 0.2 
Depreciation costs (taking into account 30 years of operation) 2.0 
Other costs 0.2 
Profit before income tax 74 
Income tax 1.5 
Net profit 5.9 


The data on annual expenses presented in Table 5 are based on the current regulatory legal acts of the Republic of 
Armenia. According to these documents, the value added tax and income tax is 20 %’, and the property tax is 0.6 % (of 
the residual value). Based on the data on the construction and current operating expenses of the SHPP, Table 6 shows 


the key indicators of financial productivity. 


’ Hayastani Hanrapetut’yan Harkayin orensgirk’, HO-165-N, yndunvats 2016 t’vakani hoktemberi 4-in (In Armenian). 
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Table 6 


Financial performance indicators 


Indicator Unit rev. Value 
Internal rate of return % 10.4 
Payback period years 9.0 


Discussion and Conclusions. Due to the mountainous terrain in the greater territory of the Republic of Armenia, 
local resistance is often created to adjust the pressure in the water supply network. As a rule, we are talking about semi- 
open valves and the installation of pressure control equipment. Equipping the network nodes with specially designed 
turbines will enable converting artificially created local resistance and extinguished mechanical energy into electrical 
energy. It should be noted that more than 300 semi-open valves and pressure control equipment were installed in the 
Yerevan water supply system alone. Many of them require replacement. The solutions proposed in this article can be 


applied to the areas under repair. 
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Abstract 

Introduction. The movement of the carrier and external factors (the effects of the atmosphere, temperature and 
pressure) degrade significantly the image quality of the servo optoelectronic systems (OES) and the positioning 
accuracy of the emitting OES. The issues of image quality improvement and the probability of keeping the image of the 
observation object (OO) on the optical axis of the servo EOS are considered. 

Materials and Methods. The development of an automatic control system for an optoelectronic device (ACS OED) 
involved solving a multi-criteria optimization problem taking into account a number of conflicting technical-and- 
economic (TE) requirements. The determination of tolerated dynamic errors (TDE) of image stabilization was a key 
issue in the development of on-board optoelectronic devices (OOED). Lagrange equations of the second kind and the 
mixed Gilbert method made it possible to obtain a mathematical model of the CO OED. Then, the decomposition of a 
two-link ACS with nonlinear cross-couplings in the CO was performed. A functional diagram of the image formation 
model of the OOED was presented. The parameters of the matrix photodetector and the requirements for the dynamic 
error of the ACS OED, taking into account the permissible MTF of the OED, were listed. The functions of transferring 
modulation, as well as linear, harmonic and vibrational shift of the image corresponding to the permissible and achieved 
TDE were visualized. Logarithmic frequency characteristics were created in the Mathcad environment. The two-link 
control system of the OED with the specified parameters of the CO for the considered movement was presented as two 
independent azimuth and elevation control channels. 

Results. The processes of control of the on-board optoelectronic system in the stabilization and tracking modes were 
described. To study the dynamics of spatial control of the OOEP in accordance with the ACS methodology, a computer 
simulation model (CSM) of the digital automatic control systems (DACS) of the OED was developed. It was 
implemented in the Matlab environment and consisted of CSM CO, drives, proportional-integral-derivative (PID) 
controllers taking into account non-linearities, a central computing device (CCD), a guidance software device, a CSM- 
carrier that implemented the equations of motion. Harmonic vibrations of the carrier were described. The errors of 
tracking and stabilization in the tracking mode with an additional control action introduced in the form of a constant 
speed were determined. The dynamics of spatial control of the OOED was investigated. A computer simulation model 
of a digital automatic control system of an optoelectronic device, the results of modeling the DACS OED without 
considering the board movement, and the processes of OED control subject to movement were visualized. 

Discussion and Conclusions. The stabilization accuracy was calculated for the studied cases. It was established that the 
stabilization tens of times exceeded the previously stated indicators, and it tens of times reduced the requirements for 


the convergence of the laser beam and the laser radiation power when developing the optical path of the product in 
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question. The proposed CSM can be used in the development of the on-board optoelectronic systems. In this case, the 
application of the presented methodology and CSM will help to reduce labor costs and minimize errors. 

Keywords: automatic control system, modulation-transfer function, servo system, optimization, optoelectronic device. 
Acknowledgements: the authors would like to thank the management of the NTI Competence Center program 
“Technologies in Robotics and Mechatronics Components” for their assistance in preparing this article. 
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Introduction. When designing modern on-board optoelectronic devices (OED) and complexes, methods of 
computer modeling of optoelectronic systems (OES) are widely used. Such authors as Yu. G. Yakushenkov, 
V. V. Tarasov, I. P. Torshina, V. P. Ivanov, V. A. Baloev, V. A. Ovsyannikov, V. L. Filippov, worked on these tasks. 
Computer modeling enables to solve the problems of rational choice of the structure, parameters, and element base of 
the OES, providing the required performance indicators under specified constraints without expensive field studies and 
tests. 

The movement of the carrier and external factors (temperature, pressure) degrade significantly the image quality of 
the servo OES and the positioning accuracy of the emitting OES. The probability of holding the image of the object of 
observation (OO) on the optical axis of the tracking OES is reduced. For emitting OES, the positioning accuracy and the 
probability of performing observation tasks are reduced. To hold the image in vision or in the irradiation sector, it is 
required to capture a fairly large solid viewing angle. Hence, it is required to increase the dimensions of optical systems 
and power electronics capacities. To avoid this, controlled radiating OES is used. The above-mentioned features of the 
OED limit the development of the ACS OED. This problem was considered by V. A. Strezhnev, V. M. Matrosov, 
A. S. Zemlyakov, N.N.Malivanov, E.I.Somov, A.I. Malikov, V.A.Krenev, A.I. Karpov, D. A. Molin, 
A. V. Mikhalitsyn. When designing controlled OED operating in guidance and tracking modes, the following fact was 
established: the guidance time and the accuracy of stabilization of the optical axis, as well as the dynamics of the OED 
subsystems affect significantly the image quality. 

Thus, the task of minimizing the time and increasing the accuracy of pointing and holding the image of the OO in 
vision of the controlled OED is a challenge. In this regard, the following is required: 

— development of mathematical models of power lines as control objects; 

— synthesis of control algorithms; 

— creation of computer simulation models (CSM); 

— study on control systems that take into account the dynamics of the movement of a controlled OED; 

— determination of parameters affecting the dynamic properties and image quality. 

The work aimed at improving the accuracy characteristics and image quality of OED operating in the guidance, 
stabilization and tracking modes, due to the rational choice of their parameters during synthesis and modeling. 

Materials and Methods. The development of ACS OED started with solving a multi-criteria optimization problem 
that took into account conflicting technical-and-economic (TE) requirements. In the development of study [1], a 
modified method of designing a self-propelled ACS OED was proposed. For the quality criterion of the ACS, we took a 
set of dynamic characteristics of control channels that met the conditions: 

Toon (Vy) = Teen (Vy), 8(@) < 64°" (w), Aa, < Aa", Ad, < Ade", (1) 
M < M*°" |A@| = Ag?”, |AL| => AL*°". 

Here, T3" (Vy), Toon(Vu.) — permissible and implemented modulation-transfer function (MTF) of the OED at the 
Nyquist frequency, 5(@); 5°" — amplitude frequency response of cross-couplings and its acceptable decomposition value; 
Aol", Adio", Aa,, Ad, — permissible and steady-state values of the dynamic errors of the ACS in angle and angular 
velocity under the action of disturbances under conditions close to the actual operation of the ACS; k= 1,2, ... — 
number of the control channel providing image quality; Myo, = 1.05—1.25 — oscillation index, Ag; AL — stability 
margins in phase and amplitude. 

The key issue in the development of the OED is the determination of tolerated dynamic errors (TDE) of image 
stabilization. 


The image formation scheme of the on-board optoelectronic device (OOED) is shown in Figure 1. 
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or Carrier model 


ACS model 
ISS model 
. device 
OO model MPD model ACC model model 


Fig. 1. Functional diagram of the OOED imaging model: OO — object of observation; OS — optical system; MPD — matrix 
photodetector; ACD — amplifying-converting device; ISS — image stabilization system 


The MTF OOED should meet the condition that provides acceptable image quality [2]: 
Tosc (v) = Tay (V)To6 W)T on W)Tyny (v)Tccu (v) > To3c (v), (2) 


(mi) | 


-1 
om = . Yu 
To5c() = 2 (sinc (w is) exp - EGa 
Here, To3c(v) — MTF OED; N — spatial frequency; T35¢(v) — admissible MTF OOED; T,,,(v) — atmospheric MTF; 
Too(v) — lens MTF; Ty,,(v) —photodetector MTF; Ty,y(v) — optical information conversion function MTF [1]; 


Tccu(v) — MTF of the image shift (dynamic error of the ISS), depending on the type of dynamic displacement of the 
image: linear (JI) — x(t) = Vt, harmonic (() — x(t) = agsin(t), random (R), and defocusing. 
The permissible MTF of the image stabilization system includes tracking, stabilization, vibration protection (VPS) 


and automatic focusing (AFS) systems [3]: 


Ta N 
Tocu(N) = Tr (N)T-(N)Tg(N) Ta (N) = TEEN) _ oe 
Ty (v) = sinc(nAay,v),T-(v) = Jo(2tAdopv), Tg(v) = exp[—2(nAagv)?], as 
2),(A he ie 
To(v) = “Ate ay = 2810 (~*) (: _ =) 


Tor(N) = Tar(N )To6(N) Ton (NT yny CN). 
Here, T;(N) — image linear shift error; T-(V) — sinusoidal error amplitude; Tg,(N) — average value of vibration 
amplitude error; J, — zero-order Bessel function of the first kind; J, (Ay) — Bessel function of the first kind of the first 
order; Ay — focus error (mm); 6 — average value of wave aberration in fractions of a wavelength; 4,, — average 


wavelength of the spectral range. 
Let us assume that each subsystem of the ACS must make the same amount of change in the image quality of 


the OED and T,,(N) = const during observation. Then, the definition of acceptable MTF R, VPS, AFS, will be 


simplified (3): 
OT yV 
TP") = * Fac G = JLT,B, ). 


We expand functions T;(N), T-CN), Tg(NV), Tp CN) (3) into a series and obtain expressions that determine the TDE: 
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aaa 0,824 j1—7A°" Gn) _ 0,335 1-72") 
Ag™ < — +, Agi™ < = —1___ (5) 


, 
Vu Vu 


0.5 in[ 73 (va))-+ 
Ag S 3,018.71 —T3°"(v,), Aan” S _—— 
Here, v,, — Nyquist frequency. 
The limiting spatial frequency that the OED should resolve during the observation process is determined by the 


Johnson criterion: 


NyL 
va Tag? V3 =05vii (6) 
where v,,— angular (bar/rad) limiting spatial frequency; Ny — number of resolution elements (Johnson 


numbers); L — distance to OO; Ay, — critical size of OO. 

When decomposing, all members of the series are discarded, except the first ones, therefore, the tolerances obtained 
should be revised by gradient descent. 

The above reasoning allowed us to determine the accuracy of image stabilization of the onboard optoelectronic 
system (OOES) for tracking and stabilization modes. A methodology and a program for calculating MTF and tolerances 
were developed (Patent RU2021660340’). 

Let us denote requirements for the dynamic error of the ACS OED, taking into account the permissible MTF OED. 
Initial data for calculation: 

— diameter of the entrance pupil — D = 90 mm, 

— OS focal length — f*= 100 mm, 

— wavelength — i= 4 um, 

— distance to the object of observation — L = 10 km, 

— meteorological visibility — SM = 15 km, 

— flight altitude — H = 5 km, 

— number of MPD elements — 256X256, 

— frame rate — F = 400 Hz, 

— probability of correct detection — P 5, = 0.8, 

— false alarm probability — P,, = 10°, 

— signal-to-noise ratio of the amplifier — optimal filter — m = 1.105, 

— angular size of the object of observation — y; = 5-10 ° rad, 

— angular distance between two adjacent targets — y, = 3.34- 10 °rad. 

MPD parameters: 

— d=0. 0.25a — thickness of the photosensitive layer, 

— €, = 4:10 ° — charge transfer inefficiency, 

— Mc = 2 — number of samples per 1 element, 

— k, = 0.6 — fill factor, 

—~1=5.4-10 °s — reading time, 

— Tap = 5-10 '°s — time constant of the transducer. 

Figure 2 shows the MTF of linear, harmonic and vibrational shift of the image corresponding to the permissible and 
obtained DDP of the ACS OOED. 


' Program for calculating the modulation transfer function of an optoelectronic device based on UAV: RF patent 2021660340, 2021. Innopolis 
University. (In Russ.) 
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Fig. 2. Modulation transfer functions T°"(v) 


In accordance with (3)—(5) and the initial data for the calculation, we determined: v,, = 600 rad' of the OOED 
observation channel in the infrared region with a detection probability of P =0.8. For MTF OOED and TDE ACS: 
Ady = 3.5 ang. min, Aor = 3.1 ang. min, Aog = 3.0 ang. min. 

On the basis of the Lagrange equations of the second kind and the mixed Gilbert method, a mathematical model of 
the SO OED [4] of guidance and tracking was obtained in the form of three information channels driven by DBM 
electric motors along the gimbal axles. We presented it in a matrix form’: 

A(q)q + B(t,q)q + F(@.q) +W,q) + 0,9) = May — Myp. (7) 
Here, q = (a 8)", a, 8 —CO rotation angles in azimuth and elevation.. 


_ (Bi + B(B) —D(B) - 0 bag(t, a, B) 
aca) = ( —D(B) Gy ).20.0= (9 cerca) 0 ) 


. _ (—24BF (a, B) — B7E(B) _ (Walt, q) 
ae ( PF (B) } wed = (vce) 
r(q) ~ ~ _ Mapa _ Myp1sign(&) 
Ota) [Gp] OAR Ma = (42s) Mo = (mn) 
bog (t, a, B) = {(2A2(B) — Cz) cos a — 2E(B) sin A} Ox, (t) - 2F (B) Wy, (t) 
— {(242(B) — Cz) sin a + 2E(B) cos o}z,(t), 
Walt, q) = {-@x, (CF (a B) + @y, (t)(By + B(B)) — wz, (t)D(a, B)} + {E (a, B)NZ, (t) — E(a, B)w?, (t) + 
+ D(a, B)a(t)w, (t) + 24s (a, Box, (ez, (t) — F(a, Bay, (Ho2, (O}, 
we(t,q) = {we, (t)[E(B) cos (a) — Cz sin(a)] — &y, (t)D(B) + wz, (CE (B) sin(@) + C, cos(a)]} - 


i 
= 5 (Fea + cos(2a)) + D(B) sin(2a)]w%,,(t) — 2F (B)w5, (0) + [F(B)(1 — cos(2a)) — D(B) sin(2a)]wz, (t) + 
+2[24,(B) cos (a) — E(B) sin(a) Joy, (t)wy, (t) — 2[F(B) sin(2a) — D(B) cos(2a)]a,, (t)z, (t) — 

—2[E (B) cos(a) + 2A2(B) sin(a@)ay, (tmz, (t)}, 

—X¢, Sin(a) + Zc, cos(a) . 
Te (@) =m, 0 + 


—X¢, COS(M) — Zc, sin(a) 


—(X¢, cos(B) — yc, sin(B)) sin(a) + Zc, cos(a) " 


+Mz 0 
—(X¢, cos(B) — yc, sin(B)) cos(a) — Zc, sin(a) 


*Burdinov KA, et al. Matematicheskaya model' i sintez sistemy avtomaticheskogo upravleniya bortovogo optiko-ehlektronnogo pribora. In: Proc. XII 
All-Russian Congress on Fundamental Problems of Theoretical and Applied Mechanics, Ufa, 2019. P. 190-192. (In Russ.) 
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(-x¢, sin(B) — Yc, cos(B)) cos(a) " 
rg(q) = mz X¢, cos(B) — Yc, sin(B) 
(xc, sin(B) + Yc, cos(B)) sin(a) 
The next stage of development was the decomposition of doubly connected ACS [5] with nonlinear cross couplings 


in the CO, which is the studied OED. The linearized equations of motion with respect to the trajectory have the form: 
@ Xzy (t) = Oot € 0 + 27; B X,, (t), = Bot € 0 + 11/2, 
(Ayo + 44, )AG + (Bry + By )AG + Cy, A + a,2AB + by AB + Cy2AB = Ky Ary, 


Az, AG + by, AG + Az ,AB + bygAB + C22AB = Ky2Au2, (8) 
Au, = R,(p)Ag,, Ae, = Aaj, — Aa*; Aug = R2(p)Ae2, Ae2 = AB, — AB”. 


Having constructed the variable coefficients of system (8), we determine their values at ft = ¢, at “dangerous” points 
at which the parameters change significantly or change signs in the range of angles (B = 0 + 7/2, a =0+ 2n). Then, 
reducing system (8) (for t = t, ) to a system with constant parameters and using Cramer's rule, we express it in the form 
of transfer functions (TF) at points t = te: 

Aa = Wi; (p, ty Au, — Wy2(p, ty )Au2, (9) 
AB = W22(p, ty Aug — Wai (p, ty) Au. 

The obtained TF (9) taking into account (8) for each selected moment of time (t= ty ) can be presented in the form of 

a block diagram with direct cross couplings (CC, Fig. 3). From now on, for simplicity of writing, we denote TF as 


Weikp) = Wikpate )- 


Fig. 3. Structural diagram of doubly connected ACS 


We evaluate the CC ACS by their frequency response (FR). For this purpose, we write down the TF of the open 
system, taking into account the 2nd closed control channel when Af, = 0 (Fig. 4). 


Wi (p) 


—> 


Fig. 4. Structural diagram of ACS, open on one control channel 
Let us build hodograph wig ©) = Ry(jo)W,.41 (jo) and tubes around it with radii for each moment of time ft, : 
E1x(@) = Ri (®) Wai (@) Wx (@) 5x (@). (10) 
Taking into account (10) and the Nyquist criterion, it is possible to judge the stability of the ACS (1), (2) and the 


influence of the CC on its stability in the range of scanning angles (3). 


As a result, for each (t, -th), we determine the CC time: 
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Wy12(@) Wx21(o) oar 
8 xi (@) a 8x12 (0) 6,.21(@)W,; (o) = ae Waal W,;(@), jij a 1,2). (1 1) 


Let us consider the CC ACS that satisfy the decomposition conditions. 
According to (8), (11), from the CO mass geometry data, for 11 selected points in time, we determined (ty = 0; 0.89; 
1.35; 1.58; 1.71; 2.3; 2.97; 3.8; 4.5; 4.75; 6.28 s). We built a,(2), bj(t), c,(d). Logarithmic frequency response (LFR) 


5,(@) (kK = 1.11) was created in the Mathcad environment (Fig. 5). 


9 | — is MELEE 
-100 — k=] 
Tees —k=2 
9 -150 | (ee 
3 i 
B -200 | | pee 
< -250 —k=8 
—k=9 

-300 | 

10° 10° 10? 10* Frequency, rad/s 


Fig. 5. Logarithmic frequency response 6,(@) 

Figure 5 shows that the entire LFR set does not exceed —50 dB. From the analysis of LFR 5,(@), we have the results 
presented below. 

YoeQ, = (10 “+10*) rad/s. 

The following was determined: max 65,(@) = 0.0136 (-37.3dB)<8"=0.0335 (-29.5 dB), 
AL = M*°"/ (M" + 1) = 0.535 (5.43 dB) at Miop= May= 1.1; WM" = 1.15. 

Hence, a doubly connected control system (8) of the OED with the given parameters of the CO for the motion in 
question can be presented as two independent control channels in azimuth and elevation angle. 

Research Results. Figures 6 and 7 show the OOED control processes in stabilization and tracking modes. The scale 
of the misalignment graphs has been increased by 600 times. The graphs of the program control and the output values 


visually coincide. 
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Fig. 6. OOED stabilization processes in azimuth 
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Fig. 7. OOED tracking processes in azimuth 


We describe the harmonic oscillations of the carrier: 4; = A.= 2 (amplitudes), 7; = 7,= 1 s (oscillation periods), and 
A,=A,= 12 T, = T,=6.28 s. We indicate the load moments: M,,; = 2.2 Nm, M,,; =0.3 Nm. With such harmonic 
oscillations and load moments, the error of stabilization of the sighting axis does not exceed: 

— in azimuth — 1.2 ang. min, 

— in elevation angle — 1.6 ang. min. 

In tracking mode, with an additional control action introduced in the form of a constant speed of 12 desiy the error of 
tracking and stabilization does not exceed 1.7 ang. min in azimuth and 1.3 ang. min in elevation angle. 

To study the dynamics of the spatial control of the OOED in accordance with the ACS technique, the CSM of the 
digital automatic control system (DACS) of a specific OED was developed (Fig. 8) [1]. 


180/3.14°60, 


beta Regulator 


OUTPUT 


, [Ch 180/3.14] { 
Software 


device 


Au 
Lie}. . ges ‘ 


Carrier Object 


Fig. 8. Computer simulation model of a digital system for automatic 
control of an optoelectronic device 


The solution is implemented in the Matlab Simulink environment. Its elements: 

— CSM CO (Object) obtained according to (7); 

— CSM of drives (Drive); 

— CSM of PID-regulators (Regulator) taking into account non-linearities and DCD; 

— CSM of software guidance device; 

— Carrier CSM that implements the equations of carrier motion. 

Figures 6 and 7 show the results of modeling the DACS OED without taking into account the movement of the 
board. In this case: 

— guidance error does not exceed 18 ang. min in azimuth, 10 ang. min in elevation angle; 

— tracking error — 0.6 ang. min in azimuth, 0.6 ang. min in elevation angle; 

— stabilization error — 0.45 ang. min in azimuth, 0.6 ang. min in elevation angle. 

Figure 9 shows the control processes of the OED, taking into account its movement. The azimuth guidance error is 
2 ang. min, the elevation angle error is 5 ang. min. The error of tracking and stabilization is 0.8 ang. min in azimuth, 


1.3 ang. min in elevation angle. 
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Fig. 9. Simulation results of the system in two-channel guidance mode (degrees): 
a — system misalignment; b — transient schedules. In second case, schedules of program control (alpha target, beta target) and 
output value (alpha, beta) visually coincide 


Discussion and Conclusions. The results obtained are presented below. 

1. A technique for evaluating tolerances for the accuracy of image stabilization of the OOED under conditions of 
controlled and temperature influences was proposed. (Patent RU2021660340). Analytical estimates of the TDE image 
stabilization according to (4), (5) for tracking and stabilization modes based on the MTF OED were presented. The 
requirements for the ACS were obtained, providing a solution to the problem of detecting OO in the infrared region. 

2. The decomposition acceptance criterion was given, the method of decomposition of a doubly connected ACS 
OEP with the use of computer technologies was described. The requirements for the design parameters of the CO of a 
real device were defined. That allowed synthesizing the ACS OEP as two independent control channels. 

3. The CSM ACS OED was proposed to study the dynamics of isolated control channels and spatial control. In this 
case, the nonlinearity and non-stationarity of a specific CO were taken into account. 

4. A mechanical model of the OED was developed, the equations of the dynamics of spatial motion (7) of the CO 
were obtained. 

5. The parameters of the PID-controllers were synthesized, providing stable guidance and stabilization of the OED 
under all specified modes. Q-factor in speed: K,, = (3700 + 4000)s71, K,, = (3000 + 3300) s~?. Corrective links: 
T11 = Tez = 0.058, Ty21 = Tya27 = 5 X 10-5 s. This provided more accurate operation of the flight planning system 
and trajectory construction [6-8]. 

6. It was shown that with the help of CSM, it was possible to obtain the required time and error of guidance, 
stabilization, and tracking. In this case, the guidance speed was 0.6 s. The corresponding indicator of foreign 
analogues — | s°. 

7. In [9], a laser counteraction system was considered. Its technical characteristics allowed performing the task at 
W= 1000 W/steradian and laser power P=3.1 W. It was shown in[10] that the required radiation power was 
W = 2700 W/steradian at a laser power of P = 200 MW. 


> Directed Infrared Countermeasure (DIRCM). EMSOPEDIA. emsopedia.org. URL: https://www.emsopedia.org/entries/directed-infrared- 
countermeasure-dircm/ (accessed: 01.05.2022). 
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We calculated the required stabilization accuracy for each of the presented cases using formula 
2x180x60 


A arccos (1 = —) [ang. min] and made Table 1. 
Table 1 


Comparison of required stabilization accuracy 


Radiation power (W), | Laser power (P), | Required stabilization accuracy (A), 
W/sterad WwW ang.min 
Installation [11] 2,700 0.2 33.4 
Installation [10] 1,000 3.1 216.1 
Product 1.3 


Thus, the accuracy of stabilization of the proposed system was ten times higher than that required in papers [9, 10]. 
And this reduces the requirements for the convergence of the laser beam and the laser radiation power tenfold when 
developing the optical path of the product in question. 


List of abbreviations and conventions 
OES — optoelectronic system, 

OED — optoelectronic device, 

OO — object of observation, 

CSM — computer simulation model, 
ACS — automatic control system, 
MTF — modulation transfer function, 
TDE — tolerated dynamic error, 

ISS — image stabilization system, 
MPD — matrix photodetector. 
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Abstract 
Introduction. Incipient fatigue damage in the metal superstructures of bridges creates certain threats to the safety of 


operation. Various methods of non-destructive testing are used for their timely detection and diagnosis. A modern and 
popular on-the-day solution is the method of infrared (IR) thermography. Due to the specifics of the operation of IR 
cameras, additional processing of recordings received from these cameras is required to obtain an accurate result. This 
work aims at presenting a method for processing thermofilms and describing the possibilities of its application under 
real conditions. 

Materials and Methods. A method for processing thermographic films was described. It provided detecting temperature 
anomalies using only information from the camera. The results of its application on the elements of existing metal 
bridge spans are presented. 

Results. It is shown that there are temperature anomalies for existing defects. This means that the defects continue to 
develop, which was confirmed by subsequent observations of their condition. In addition, a case of temperature 
anomaly in the defect-free external region was identified. This might be a sign of an incipient defect that could not be 
diagnosed by other methods. If the presence of this defect is confirmed during repeated examinations, it will be possible 
to diagnose hidden defects that have not yet come to the surface, and/or detect potentially collapsing places. 

Discussion and Conclusions. The IR thermography performance as a method of non-contact non-destructive testing is 


shown, as well as its operability on real objects under random load. 
Keywords: IR thermography, nondestructive testing, fatigue cracks, metal bridges, structural defects, IR camera. 
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Introduction. Fatigue cracks in the metal structures of bridges, formed under the influence of transport loads, wind 
loads, etc., are a serious challenge. Although these cracks in themselves only affect the durability of the structure, under 
certain conditions, they can give rise to the development of brittle cracks, which directly affects the safe use of 


structures [1]. 


© Solovyev A. L., Royak M. E., 2022 
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One of the modern methods of non-destructive testing, which is gaining popularity, is IR thermography [2-12]. The 
method principle is based on the fact that any body that is subjected to mechanical action releases an amount of heat 
proportional to the intensity of the action. Therefore, when examining a section of the structure containing defects, you 
can immediately get information about the distribution of stresses in this section. A different temperature distribution 
will be an indicator of ongoing fatigue damage [13, 14]. 

The major challenge when working with IR cameras (thermal imagers) is to obtain results with acceptable accuracy. 
Due to the structural features of the detectors, the highest accuracy is achieved when using cold sensor cameras. 
However, the high cost of such devices is a serious obstacle to their use. The introduction of cameras using uncooled 
arrays of microbolometers improved the situation, as these cameras turned out to be noticeably cheaper. The 
disadvantage of such thermal imagers is the lower accuracy compared to cameras with cooled matrices. Therefore, to 
get stable results, additional processing of the obtained thermograms is required. 

This work aimed at evaluating the effectiveness of self-reference lock-in thermography using a signal obtained 
directly from thermal imaging’. 

Materials and Methods. The IR thermography method is based on the analysis of the temperature distribution in 
the material. When tensile or compressive stresses are applied to a material, its temperature changes in proportion to the 
magnitude of these stresses in accordance with the following formula”: 

AT =-—*~_T Ao: (1) 
pC, 
where o, — thermal-expansion coefficient; p — matter density; C, — specific heat capacity at constant pressure; T— absolute 
temperature, K; Ao — change in principal stresses, MPa. 

When adiabatic conditions are maintained, as well as when elastic deformations develop only, the coupling (1) is 
linear and reversible, in which a change in temperature follows a change in stress. As soon as damage begins to form in 
the material, it ceases to be elastic and develops plastic deformations. The energy of a system subject to such 
deformations is intensively converted into thermal energy, thus making the main contribution to the thermal effect. 
Hence, the thermal picture becomes a reflection of the internal state of the material (the degree of destruction) and, 
consequently, an indicator of the stress level. 

The relative amplitude of an IR signal can be determined using the least squares method®. Let us assume that load F, 
whose impact in a certain area can be measured and denoted as f(t), acts on the structure. We will consider this signal as 
the reference one. Suppose there is also an IR video recording of this construction in the same time interval, which can 
be presented as a sequence of frames (a frame is a matrix of values M x K in size) in time, and we denote values (i, j) of 


the element of this matrix (pixel) at time t where 1<i<M, 1<j<K, as y(t). We denote the approximation 


function for this element as: 


Y,(t)=a, +b, f(t). (2) 


where a, — bias; b, — coefficient of influence of the reference signal f. 


*Sakagami T, Nishimura T, Kubo S, et al. Development of a Self-reference Lock-in Thermography for Remote Nondestructive Testing of Fatigue 
Crack (1st Report, Fundamental Study Using Welded Steel Samples). Transactions of the Japan Society of Mechanical Engineers. Series A. 
2006;72:1860-1867. https://doi.org/10.1299/kikaia.72.1860 

Thomson W. (Lord Kelvin). On the dynamical theory of heat. Trans. R. Soc. Edinburgh. 1853;20:261-283. 

> Lesniak JR, Boyce BR, Howenwater G. Thermoelastic Measurement Under Random Loading. In: Proc. SEM Spring Conf., 1998. P. 504-507. 
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Using function (2), we approximate signals y(t) from the thermal imaging film. Let us show how to find these 
coefficients for some y;,(t), the reasoning for the rest will be similar. We will minimize the deviation of our 


approximation from the signal received from the camera: 


N 


A =3(y,(n)-¥,(n)f > min. @) 


n=l 


where N — number of frames in the recording; n — frame number corresponding to the time. From here, b, can be 


found as follows: 


- ¥ 9u(n) 
YI) Youle) WL vy lod elo)- DAIS v4 (0) 5 
, N > f(n) nS r°eo)-(S 10) 


Through repeating (3) and (4) for all (i, j), we get matrix B= 1b, j of the same size as the original shooting frame. 


ij 
Thus, we obtain an approximation of each signal y,(t) of the IR survey in time by means of reference signal f(t). The 
values of matrix B show the relative intensity of the temperature change in a certain area compared to the intensity of 
the reference signal change. Reference signal f(t) may not necessarily be received from a third-party system, e.g., froma 
tensometer. The described approach, called self-reference lock-in thermography, can also be applied with a reference 
signal obtained from the same IR recording* [15, 16]. 

Research Results. The proposed algorithm was implemented in the form of Python scripts. With the help of these 
scripts, thermographic films recorded on the metal spans of existing automobile bridges, in whose structures fatigue 
cracks had been previously diagnosed, were processed. The performance of the algorithm was verified on the bench 
tests simulating the behavior of real structures with known information about the defect. The survey was carried out 
using an infrared camera with uncooled microbolometer Fluke Ti 400 having a thermal sensitivity of less than 0.05°C 
and a shooting frequency of 9 Hz. The recordings were made at the time of the impact of the automobile load on the 
bridge. As reference signals for each record, a defect-free zone of 15x15 pixels was used from the same structural 
element near the existing and/or proposed destruction zone. Processing results were visually presented in the form of 
images obtained by means of matrix B values. Each value of the matrix corresponded to one pixel in the image, the 
number was converted to color using the palette function defined in the matplotlib package. 

Figure | shows a fragment of a superstructure steel beam. A crack was found at the junction of the horizontal sheet 
with the wall (Fig. 1 c). Figure 1a shows one of the frames of IR film of this crack. Figure 1 b presents the result of 
processing, which clearly shows a spot, indicating heating occurring in the designated area. This zone corresponds to 
the area around the crack tip, which indicates the continued development of the defect. This was also confirmed by 


repeated surveys, which established an increase in cracks of more than 30% over 4 years. 


“ Galietti U, Modugno D, Spagnolo L. A novel signal processing method for TSA applications Measurement Science and Technology. 2005;16:2251. 
http://dx.doi.org/10.1088/0957-0233/16/11/017 
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c) 
Fig. 1. Infrared results: 
a — film frame; b — image built on the values of matrix B; c — photograph of the crack 


Figures 2 and 3 show cracks that originated in the weld, then bifurcated and went into the stiffener and decking. The 
survey was carried out under the bridge tests. Strain gauges were installed in the same places to monitor the stresses, the 


results of which confirmed the development of these cracks. 


Fig. 2. Diagnosed fatigue crack: 
a — image built on the values of matrix B; b — photograph 


Reference marks 
for image alignment 
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c) 
Fig. 3. Diagnosed fatigue crack: 


a — film frame; b — image built on the values of matrix B; c— photograph of a crack 
Figure 4 shows the combined readings of tensometers and temperature for the same section of the structure. 
Temperature indicators were obtained through averaging in the survey area and subsequent smoothing over time. The 
strain gauge data were also smoothed. As can be seen, the temperature change completely repeats the voltage change up 


to the sensitivity of the device. 


Normalized values, ea 


——— Voltage change 


0.1 Temperature change | 
Ol L tL L _ 
5 10 15 20‘ Time, s 


Fig. 4. Voltage changes graphs combined with temperature changes (temperature changes are taken with the opposite sign). 
Values along the y-axis are normalized with respect to unity 


Figure 5 shows a structural element that had no visible damage, which was also not unambiguously determined by 
other methods. However, as a result of processing the thermogram package, it was found that there was a thermal 
anomaly at the junction of two elements. This indicates that the material was probably self-heating at this point, and a 
defect was emerging under the surface, which would soon have to come to the surface. Thus, this section of the 


structure requires additional attention under future surveys. 
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a 


a) b) 
Fig. 5. Infrared results: a — weld seam without visible fatigue damage; b — image of the weld toe, 
built from the values of matrix B 


Discussion and Conclusions. The work carried out on the superstructures showed the effectiveness of self- 
reference lock-in thermography as a method of contactless diagnostics. The method performance on real objects under 
random load was shown. There was a correspondence of the stress change in the structure with the change in the 
recorded temperature. We carried out the diagnostics of some known defects, on which the continued development of 
damage was noted, which was subsequently confirmed by repeated examinations of these structural elements. In 
addition, a case of self-heating in an externally defect-free zone was identified, which indicated a probable process of 
the origin of the defect, which at the time of shooting could not be unambiguously diagnosed by other methods. This 
proves that this site requires additional monitoring. If a defect is detected at this site in the course of subsequent 
observations, we can talk about the possibility of using IR thermography to identify hidden defects and predict their 
appearance. 
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Abstract 
Introduction. The article presents an overview of modern neural network models for natural language processing. 


Research into natural language processing is of interest as the need to process large amounts of audio and text 
information accumulated in recent decades has increased. The most discussed in foreign literature are the features of the 
processing of spoken language. The aim of the work is to present modern models of neural networks in the field of oral 
speech processing. 

Materials and Methods. Applied research on understanding spoken language is an important and far-reaching topic in 
the natural language processing. Listening comprehension is central to practice and presents a challenge. This study 
meets a method of hearing detection based on deep learning. The article briefly outlines the substantive aspects of 
various neural networks for speech recognition, using the main terms associated with this theory. A brief description of 
the main points of the transformation of neural networks into a natural language is given. 

Results. A retrospective analysis of foreign and domestic literary sources was carried out alongside with a description 
of new methods for oral speech processing, in which neural networks were used. Information about neural networks, 
methods of speech recognition and synthesis is provided. The work includes the results of diverse experimental works 
of recent years. The article elucidates the main approaches to natural language processing and their changes over time, 
as well as the emergence of new technologies. The major problems currently existing in this area are considered. 
Discussion and Conclusions. The analysis of the main aspects of speech recognition systems has shown that there is 
currently no universal system that would be self-learning, noise-resistant, recognizing continuous speech, capable of 


working with large dictionaries and at the same time having a low error rate. 
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semantic consistency. 
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Introduction. This article provides an overview of the main language model based on the neural network for 
Natural Language Processing that helps computers communicate with people in their native language and scale other 
language tasks. Modern machine learning technologies allow computers to read text, hear speech, interpret it, 
measure moods, and determine which parts of speech are important. This technology is called Natural Language 
Processing (NLP), it is based on many disciplines, including computational linguistics. NLP is increasingly being used 
in interactivity and productivity applications, such as creating spoken dialogue systems and speech-to-speech engines, 
searching social networks for health or financial information, detecting moods and emotions towards products and 
services, etc. 

The relevance of NLP is primarily associated with the need to process large amounts of audio and text information 
accumulated by mankind over the past decade. Currently, most modern devices are endowed with a voice control 
function, and various kinds of digital assistants are becoming widespread. Now, the speech recognition function is 
available in almost any gadget, it allows us to interact through voice applications, facilitating and simplifying a person's 
life. There are a fairly large number of commercial speech recognition systems, among the most famous there are 
Google, Yandex, Siri. The quality of speech recognition in such systems is at a fairly high level, but they are not 
without a number of shortcomings. Unfortunately, despite the amazing development of computer technology, the 
current problem of equipping a computer with a full-fledged, natural human voice interface is still far from over. 

Materials and Methods. NLP technology is rapidly advancing due to the increased interest in the field of machine 
learning, as well as the availability of big data, powerful computing, and improved algorithms. However, NLP is not a 
new science. Attempts to teach computers to communicate with people through a natural voice interface have been 
made since the early days of computer technology. NLP was born with the advent of the first computers from the idea 
of how good it would be to use these machines to solve various useful tasks related to natural language, e.g., these 
programs were intended for people who, due to physiological characteristics, could not type text manually [1]. 

The first task that the first computers solved in the early stages of the formation of NLP was the task of machine 
translation, i.e., automatic translation of text from one language to another using a computer. This problem was 
successfully solved and started to be applied in the mid-1950s, in the past century, for the “Russian-English” pair [2]. 

The second task of machine learning was to create conversational systems, the programs to conduct a dialogue with 
a person in natural language. Many systems created at that time were imperfect due to a number of difficulties in speech 
recognition that can have a significant impact on the quality of the result. The difference in the voices of speaking 
people, the inconsistency of colloquial speech, the phonogram of the same words can vary greatly depending on a 
number of factors: pronunciation speed, regional dialect of the language, foreign accent, social class, and even the 
gender of a person [3]. 

The third task was to create a question-and-answer system. There was a need for programs that would answer 
exactly the human question. At that stage, such a question was in the form of a natural language text. Thus, the problem 
of scaling the recognition system has always been a significant obstacle. In the course of many years of research, it has 
been found that it is required to involve not only programmers, but also experts in linguistics, radio engineers, 
mathematicians, biologists, and even psychologists in solving the problem. 

At different times, various mathematical, statistical, logical, stochastic approaches were used in natural language 
processing, such as Dynamic Time Warping, Bayesian discrimination, Hidden Markov Model, formal grammars, and 
probabilistic approaches. At the present stage of natural language processing, machine learning methods are 
widespread, in particular, neural networks. Currently, in modern linguistic research, at the first stage, texts are selected 
that are planned to be analyzed, and a corpus of texts is created. Next step, the collected material is transferred to an 


expert linguist. He prescribes the rules, compiles dictionaries, marks up texts for the identification of target structures in 
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texts for the further task solution. Another method is also used, in which an expert linguist marks the text into target 
structures or categorizes texts in certain classes, and then machine learning methods automatically derive some rules or 
models for further solving current problems. At the end of the work, the quality of the methods is always checked. 

Philologists study semantics of the text considering meanings of polysemantic units in context, emphasizing that 
context plays a fundamental role in the word definition. Therefore, e.g., the authors discover the contextual meanings of 
polysemantic units that are not registered in lexicographic sources. In the early stages, scientists proposed to divide any 
sentence into a set of words that could be processed individually, which was much easier than processing a whole 
sentence. This approach is similar to the one used to teach a new language to children and adults. When we first start 
learning a language, we are introduced to its parts of speech. Let us consider English as an example. It has 9 main parts 
of speech: noun, verb, adjective, adverb, pronoun, article, etc. These parts of speech help to understand the function of 
each word in a sentence. However, it is not enough to know the category of a word, especially for those that may have 
more than one meaning. Specifically, the word “leaves” can be the verb “to leave” in the 3rd person singular or the 
plural form of the noun “leaf”, which should be considered from the point of language as a system of interrelated and 
interdependent units. The idea of consistency in the lexical and semantic sphere of language was first expressed by 
M. M. Pokrovsky, emphasizing that “words and their meanings do not live a separate life from each other, but are 
connected (in our soul), regardless of our consciousness, into different groups, and the basis for their grouping is 
similarity or direct opposition in their basic meaning”. Paradigmatic, syntagmatic, and epidigmatic relations among 
language units are important manifestations of the systematic and regular nature of language. The researchers note that 
words enter the syntagmatic relations based on the logical contiguity of concepts and, consequently, their compatibility 
with each other [4]. 

We need to understand that from the point of view of computer science, speech is not structured information, but a 
sequence of characters. To ensure that voice data can continue to be used, the speech recognition application translates 
it into text. The accent, individual intonations, and emotions are already being erased in the text. When data are 
translated into text, they are translated with zeros and ones. 

Therefore, computers need a basic understanding of grammar to refer to it in case of confusion. Thus, the rules for 
the structure of phrases appeared. They are a set of grammar rules by which a sentence is constructed. In English, it is 
formed with the help of a nominal and a verb group. Consider the sentence, “Kate ate the apple”. Here, “Kate” is a noun 
phrase and “ate the apple” is a verb phrase. Different sentences are formed using different structures. As the number of 
phrase structure rules increases, a parse tree can be created to classify each word in a particular sentence and arrive at its 


general meaning (Fig. 1). 


Noun words group 


The spoken language system integrating speech recognition and speech synthesis is the core technology of human- 


SENTENCE 


Verb words group 


Fig. 1. Parse Tree 


computer interaction, and oral understanding is the core of spoken language system [5]. 
Artificial neural networks, created in the form of computer models, cope successfully with the tasks of pattern 


recognition. They are trainable, they can be easily adapted (and this has already been done) to solve many 
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practical problems related to speech recognition, control of various machines and devices, event prediction, etc. 
The biggest feature of a deep neural network is training a large amount of data, and then extracting characteristic 
information. Characteristic information obtained through a network structure can give good results in the speech 
comprehension tasks. 

An advantage of a neural network in speech processing is that the perceptron can perform discriminant learning 
between speech units that represent the output classes of the perceptron. The perceptron not only learns and optimizes 
the parameters for each class on the data belonging to it, but also tries to reject the data belonging to other classes. The 
perceptron is a structure with a high degree of parallelism, which allows the use of parallel hardware. The first neural 
network models used in speech recognition systems were developed for static signals, and not for their sequences or 
signals subject to temporal variability. Later, recurrent neural networks and convolutional neural networks were 
proposed [5]. 

As language models are trained on larger and larger texts, the number of unique words (the vocabulary) increases, 
and continuous space embedding helps to alleviate the curse of dimensionality in language modeling. An alternate 
description is that a neural net approximates the language function. The neural net architecture might be fed forward or 
recurrent, and while the former is simpler, the latter is more common. Improved algorithms, powerful computers, and 
an increase in digitized data have fueled a revolution in machine learning. And new techniques in the 2010s resulted in 
“rapid improvement in tasks” including language manipulation, in particular, transformer — architecture based on a deep 
learning model. It was first introduced in 2017 [6]. Table 1 presents the main types of currently existing language 
models of neural networks. 


Table | 
Main types of neural network language models 


Language model Characteristics 


Bidirectional Encoder Representations from Transformers is a new method of pretraining 
BERT-base (2018) language. BERT is different because it is designed to read in both directions at once. Using 
this bidirectional capability, BERT is pretrained on two different, but related, NLP tasks: 
Masked Language Modeling and Next Sentence Prediction [7, 8]. 


Embeddings from Language Model is a word embedding method for representing a sequence 
of words as a corresponding sequence of vectors, but unlike BERT, the word embeddings 


ELM 
(201 . produced by the “bag-of-words” model is a simplifying representation. ELMo embeddings are 
context-sensitive, producing different representations for words that share the same spelling 
but have different meanings (homonyms) [9]. 
GPT is a Transformer-based architecture and training procedure for natural language 
GPT (2018) processing tasks. Training follows a two-stage procedure. First, a language modeling objective 


is used on the unlabeled data to learn the initial parameters of a neural network model. 
Subsequently, these parameters are adapted to a target task using the corresponding supervised 
objective [10]. 


ESPnet mainly focuses on end-to-end automatic speech recognition (ASR), and adopts widely- 
used dynamic neural network toolkits, Chainer and PyTorch, as a main deep learning engine. 
ESPnet (2018) ESPnet also follows the Kaldi ASR toolkit style for data processing, feature extraction/format, 
and recipes to provide a complete setup for speech recognition and other speech processing 
experiments [11]. 


Model uses only ID convolutions, batch normalization, ReLU, dropout, and residual 


201 
Jasper (2019) connections [12]. 


GPT-2 translates text, answers questions, summarizes passages, and generates text output on a 
level that, while sometimes indistinguishable from that of humans, can become repetitive or 
nonsensical when generating long passages [13]. 


GPT -2 (2019) 


WAV2LETTER++ It is an open-source deep learning speech recognition framework. wav2letter++ is written 


(2019) entirely in C++, and uses the ArrayFire tensor library for maximum efficiency [14]. 
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Language model Characteristics 
WAV2VEC (2019) wav2vec, is a Comvelnatal eet network that takes a audio as input and computes a 
general representation that can be input to a speech recognition system [15]. 
These are cross-lingual language models (XLMs): one unsupervised that only relies on 
monolingual data, and one supervised that leverages parallel data with a new cross-lingual 
XLM (2019) air ; : ae 
language model objective. It obtains state-of-the-art results on cross-lingual classification, 
unsupervised and supervised machine translation [16]. 
XLNet uses a generalized autoregressive retraining method that enables learning bidirectional 
XLNet contexts through maximizing the expected likelihood over all permutations of the factorization 
(2019) order and autoregressive formulation. XLNet integrates ideas from Transformer-XL, the state- 
of-the-art autoregressive model, into retraining [17]. 
This implementation is the same as BERT Model with a tiny embeddings tweak as well as a 
RoBERTa (2019) setup for ROBERTa pretrained Sa RoBERTa has the same eee as on pe 
uses a byte-level BPE as a tokenizer (same as GPT-2) and applies a different pretraining 
scheme [18]. 
ELECTRA Efficiently Learning Encoder That Accurately Classifies Token Replacements is a new pre- 
(2020) learning method that outperforms development estimation without increasing the 


computational cost [19]. 


STC system aims at multi-microphone multi-speaker speech recognition and diarization. 
The system utilizes soft-activity based on Guided Source Separation (GSS) front-end and 
STC System (2020) a combination of advanced acoustic modeling techniques, including GSS-based training 
data augmentation, multi-stride and multi-stream self-attention layers, statistics layer and 
spectral [20]. 


Unlike other models created to solve specific language problems, their API can solve “any 
GPT — 3 (2020) problems in English”. The algorithm works on the principle of autocompletion: you enter the 
beginning of the text, and the program generates the most likely continuation of it [21]. 


ALBERT incorporates two parameter reduction techniques that lift the major obstacles in 
scaling pretrained models. The first one is a factorized embedding parameterization. By 
splitting a large vocabulary embedding matrix into two small matrices, it separates the size of 


el BEE ese) the hidden layers from the size of vocabulary embedding. The second technique is cross-layer 
parameter sharing. This technique prevents the parameter from growing with the depth of the 

network [22]. 
Pretrained BERT with Whole Word Masking due to the complexity of Chinese grammar 
BERT-wwn-ext, structure and the semantic diversity, a BERT (wwm-ext) was proposed based on the whole 
2021) Chinese word masking, which mitigates the drawbacks of masking partial Word Piece tokens 


in pretrained BERT [23]. 


This is Pathways Language Model 540-billion parameter, dense decoder. Only Transformer 
PaLM (2022) model trained with the Pathways system enabled us to efficiently train a single model across 
multiple TPU v4 Pods [24]. 


As can be seen from Table 1, the first transformer models, using a bidirectional capability, allowed two different but 
related tasks of the NLP to be studied beforehand: simulating a masked language and predicting the next sentence. 
Bidirectional Encoder Representations from Transformers consist of two steps: the first step is pretraining where the 
data enter the layers of transformer, and the result of this step are vectors for words. The second step is fine tuning. The 
pretraining step consists of two steps: the masked LM and Next Sentence Prediction (NSP) [7, 8]. BERT is not without 
flaws, the most obvious one is the learning method — the neural network tries to guess each word separately, which 
means that it loses some possible connections between words during the learning process. Another one is that the neural 
network is trained on masked tokens, and then used fundamentally different tasks, more complex ones. 

Embeddings from Language Model is a deep contextualized word representation that models both complex 
characteristics of word usage (e.g., syntax and semantics), and how this usage varies across linguistic contexts (i.e., to 


model polysemy), such as “bank” in “river bank” and “bank balance”. These word vectors are learned functions of the 
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internal states of a deep bidirectional language model (biLM), which is pretrained on a large text corpus. They can be 
easily added to existing models and significantly improve the state of the art across a broad range of challenging NLP 
problems, including question answering, textual entailment, and sentiment analysis [9]. 

To alleviate the problem, suffering from the discrepancy between the pretraining and fine-tuning stage because the 
masking token [MASK] never appears on the fine-tuning stage, XLNet was proposed, which is based on Transformer- 
XL. To achieve this goal, a novel two-stream self-attention mechanism, and one to change the autoencoding language 
model into an autoregressive one, which is similar to the traditional statistical language models, were proposed [17]. 
RoBERTa, STC System, GPT models were used in quite a large number of systems. And they showed pretty good 
results. These models suggested that averaging all token representations consistently induced better sentence 
representations than using the token embedding; combining the embeddings of the bottom layer and the top layer 
outperformed the use of the top two layers; and normalizing sentence embeddings with a whitening algorithm 
consistently boosted the performance [18, 20, 21]. 

The next step, probably, will be to study the oversampling and undersampling of textual data to improve the overall 
entity recognition effect. 

Results. The analysis of the literary sources describing new methods of processing oral speech, which provides 
information about neural networks, methods for the structure and synthesis of speech, made it possible to detect 
the following: 

1. All the models presented in the review require large computing power to solve natural language processing 
problems. It is computationally more expensive due to its larger structure. 

2. None of the currently existing technologies enable solving the full range of tasks for recognizing continuous, 
defective speech. 

3. Most natural language processing models are designed to handle a wide variety of English dialects and idioms. 

Discussion and Conclusions. Voice assistants reproduce and reinforce all stereotypes algorithms. They, as a rule, 
reproduce those stereotypes that exist now in society. What does this achievement really mean? It means that the voice 
assistant is no worse (or maybe even better) than an average person at recognizing the speech of a person with a 
standard North American accent. But if an African American speaks to an assistant, then the accuracy will drop to about 
80 %. This is a huge difference. Moreover, when converting voice to text, the specifics of writing, which can be 
important for speakers, are guaranteed to be lost. 

Voice assistants do not take into account the speech and user habits of the elderly and people with special needs. 

And here, it is not even always the complexity of recognition. There is, e.g., such a condition as dysarthria — a 
feature of the functioning of the connections between the speech apparatus and the nervous system, which can cause 
difficulties in pronouncing individual sounds or, in general, in speech. 

Also, due to hardware limitations, any cartridge will result in too many model parameters and unsuccessful 


execution. The way to solve the problem of multiple cycles of dialogue requires further research. 
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